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1 Introduction 

In any experiment, external noise sources and control errors limit the accuracy of the prepa- 
ration and manipulation of quantum states. In quantum computing, these effects place an 
important fundamental limit on the size and accuracy of quantum processors. These restric- 
tions may be reduced by quantum error correction. Although very sophisticated quantum 
error-correcting codes exist which are robust against any general error, these codes require 
large-scale multipartite entanglement and are very challenging to implement in practice [1]. 
Therefore, it is of great interest to investigate schemes which reduce errors with a smaller 
resource overhead. One alternative strategy involves replacing an error-prone operation by 
a pulse sequence which is robust against the error. 

The basis of our strategy is that all noises and errors can be treated as an unwanted 
dynamic generated by an error Hamiltonian. This error Hamiltonian can arise either through 
interactions with the environment or by the misapplication of control fields. This view unifies 
the pulse sequences developed for combating unwanted interactions, e.g., dynamic decoupling 
and dynamically corrected gates [2, 3, 4], with those for overcoming systematic control errors, 
compensating composite pulse sequences [5]. In each case, the methods are limited by the 
rate at which control occurs relative to the time scale over which the error Hamiltonian 
fluctuates. However, many experiments are limited by control errors and external fields that 
vary slowly relative to the time scale of a single experimental run but vary substantially over 
the number of experiments required to obtain precise results. 

In this review paper, we examine a number of techniques for handling unwanted control 
errors including amplitude errors, timing errors, and frequency errors in the control field. We 
emphasize the common principles used to develop compensating pulse sequences and provide 
a framework in which to develop new sequences, which will be of use to both quantum 
computation and coherent atomic and molecular spectroscopy. 



2 COHERENT CONTROL OVER SPIN SYSTEMS 3 

2 Coherent control over spin systems 

The accurate control of quantum systems is an important prerequisite for many applications 
in precision spectroscopy and in quantum computation. In complex experiments, the task 
reduces to applying a desired unitary evolution using a finite set of controls, which may 
be constrained by the physical limitations of the experimental apparatus. Stimulated by 
practical utility, quantum control theory has become an active and diverse area of research 
[6, 7, 8, 9, 10]. Although originally developed using nuclear- magnetic resonance (NMR) 
formalism, compensating pulse sequences can be approached from the perspective of quantum 
control theory with unknown systematic errors in the controls [11]. Here we review several 
fundamental concepts in quantum control, and apply these ideas using NMR as an instructive 
example. Although we restrict the discussion to control in NMR spectroscopy, the following 
analysis is quite general as several other coherent systems (e.g. semiconductor quantum dots 
[12], superconducting qubits [13, 14, 15], and trapped ions [16, 17]) may be considered by 
minor modifications to the Hamiltonian. 

In practice, a desired evolution is prepared by carefully manipulating the coupling of the 
system to a control apparatus, such as a spectrometer. In the absence of relaxation, the 
coherent dynamics are governed by the quantum propagator U{t), which in non-relativistic 
quantum mechanics must satisfy an operational Schrodinger equation. 



-i(y2u,{t)HA 



U{t) = -ilJ2u,{t)HAu{t), f/(0) = l. (1) 

In this model, the unitless Hamiltonians H^ G {Hi,H2, . . . ,Hn} are modulated by real- 
valued control functions u^(t) G {ui{t) , U2{t) , . . . , u„(t)} and represent the n available degrees 
of control for a particular experimental apparatus. In analogy with linear vector spaces, we 
interpret the vector u(t) = {ui(t) , U2(t) , . . . , Unit)) as a vector function over the manifold 
of control parameters, with components Uf^{t) representing the magnitudes of the control 
Hamiltonians with units of angular frequency. It is convenient to introduce a second vector 
of control Hamiltonians H = [Hi, Hi, . . . , if„) and the short-hand notation H{t) = u{t) H. 
We omit a term which represents the portion of the total Hamiltonian which is outside of 
direct control (i.e. a drift Hamiltonian); in principle it is always possible to work in an 
interaction picture where this term is removed. Alternatively, one may assign a Hamiltonian 
Hq to represent this interaction, with the understanding that Uo{t) = 1 for all t. 

For a given control system, a natural question concerns the optimal approximation of a 
desired unitary propagation using a set of constrained control functions. Constraints may 
include limitations on the total operation length, control amplitudes or derivatives. The 
study of this question requires the solution of Eq. (1) for a particular set of controls; such 
solutions may be obtained using several methods, including the Dyson series [18], and the 
Magnus [19, 20], Fer [21], and Wilcox [22] expansions. We label particular solutions to the 
control equation over the interval ti < t < tj as U{u(t);tf,ti). If the set of all possible 
solutions to Eq. (1) is the set of all unitary gates on the Hilbert space (i.e. the solutions 
form a representation of U(n) or SU(ra), to be discussed in section 3.1) then the system is 
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Figure 1: The experimental controls available to manipulate a quantum system are mod- 
eled using a set of real- valued control functions {u^(t)} that modify a set of dimensionless 
Hamiltonians {i?^}- a) An example control function u^{t). b) A discrete approximation for 
u^{t) composed of square pulses. 

operator controllable [23, 10]. 

For some applications, it is convenient to assume that the operation time r is discretized 
into m-many time intervals, over which the control functions u^ are constant, i.e., during the 
/cth time interval Atk the Hamiltonian is time-independent and the resulting evolution oper- 
ator is Uk = U{uk] tk + Atk,tk) = exp(— iAtfc-Ufc ■ H). Figure 1 illustrates an example control 
function and a possible discretization scheme. Over each time interval, the applied unitary 
operation is a square pulse. In many experiments, the application of gates using sequences 
of square pulses is preferred for simplicity. The discretization of the control functions into 
square pulses allows one to solve Eq. (1) in a piecewise fashion. The total propagator for a 
sequence of m-many time steps is given by the time-ordered product 



U{u{t), r, 0) = t/^^™_i ...Uk...U2U, = l[Uk 



(2) 



k=l 



where the multiplication of each successive operator is understood to be taken on the left; 
this is in agreement with standard quantum mechanics conventions, where operations are 
ordered right-to- left, but at odds with some NMR literature where successive operations are 
ordered left-to-right. Frequently we will consider propagators over the entire duration of a 
sequence < t < r; as a matter of notational convenience we drop the time interval labels 
whenever there is no risk of confusion. If a pulse sequence U{u{t)) is equivalent to a target 
operation Ut in the sense that during an experiment, U{u{t)) may be substituted for Ut, 
then U{u{t)) is a composite pulse sequence [24, 25]. 

The usefulness of composite pulse sequences lies in that in many cases, one may simulate 
a target unitary transformation Ut, which may be difficult to directly implement, by instead 
implementing a sequence of simpler pulses which under some set of conditions is equivalent 
to Ut- Composite pulse sequences may be designed to have several important advantages 
over a directly applied unitary, such as improved resilience to errors. The properties of pulses 
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or pulse sequences may be considered either by the transformations produced on a particular 
initial state p, or by comparing the sequence to an ideal operation, which contains information 
on how the sequence transforms all initial states. Following Levitt [5], we assign composite 
pulses into two classes: the fully-compensating class A, and the partially-compensating class 
B. The properties of these classes are briefly reviewed. 

Class A: All composite pulse sequences in class A may be written in the form 

e^i^Uiuit)) = Ut, (3) 

where it is assumed that the individual pulses in the sequence are error- free. The global phase 
(f) is irrelevant to the dynamics since for any initial state p, both the pulse sequence and the 
target operation apply the same transformation {U {u{t)) pU\u{t)) = UxpU^)- Sequences in 
this class are suited for use in quantum computation, since the transformation is independent 
of the initial quantum state. The study of these sequences will be the primary topic of this 
article. 

Class B: Composite pulse sequences in class B transform one particular initial condition 
to a set of final conditions, which for the purposes of the experiment, are equivalent. For 
example, consider an NMR experiment on a spin I = 1/2 nucleus where the nuclear magne- 
tization, initially oriented "spin-up" (i.e. p = 1/2 + H^, where H^, Hy, H^ are the nuclear 
angular momentum operators) is transferred to the Hx-Hy plane by a pulse sequence. In 
this sense, all sequences that apply the transformations are equivalent up to a similarity 
transform exp(— i/3ifz), which applies a H^ phase to the spin. Other sequences in class B 
may satisfy U {u{t)) pU\u{t)) = UxpU]^ for a particular initial state p, but fail to satisfy Eq. 
(3). All sequences in class B are generally not well suited for use in quantum computation, 
since implementation requires specific knowledge of the initial and final states of the qubit 
register. Class B sequences are however very useful in other applications, including NMR 
[26, 24, 25], MRI [27], control over nitrogen vacancy centers [28] and in ion trapping experi- 
ments [29]. We close our discussion by noting that class B sequences may be converted into 
a fully-compensating class A sequence by a certain symmetrical construction [30] . 

2.1 Errors in quantum control 

We now consider the effects of unknown errors in the control functions. Recall that a 
pulse sequence may be specified by a set of control functions {^^(t)}, which we group into 
a control vector u{t). Suppose however, during an experiment an unknown systematic 
error deforms each of the applied controls from u^{t) to v^{t). In the presence of unknown 
errors, the perfect propagators U{u{t);tf,ti) are replaced with their imperfect counterparts 
V{u(t); tf, ti) = U{v{t); tf,ti), which may be regarded as an image of the perfect propagator 
under the deformation of the controls. 

In practice, systematic errors typically arise from a miscalibration of the experimental 
system, for example an imprecise measurement of the intensity or frequency of a controlling 
field. In these cases, it is appropriate to introduce a deterministic model for the control de- 
formation. Specifically, v{t) must be a function of the controls u{t), and for each component 
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Figure 2: Systematic errors induce deformations of the ideal control functions Uf^{t) (solid 
curves) to the imperfect controls v^(t) = u^{t) + e5u^{t) (dashed curves). Common error 
models include a) a constant unknown offset in the control, and b) an error in the amplitude 
of the control function. 



of the imperfect control vector we may write 



v,,{t) = f^[u{t);e\, 



(4) 



where the functional /^ is the error model for the control and the variable e is an unknown 
real parameter that parametrizes the magnitude of the error. This construction may be 
generalized to the case of multiple systematic errors by considering error models of the form 
f^[u(t)] ei,ej, . . . ,ek]. The model /^ is set by the physics of the problem, and is chosen to 
produce the correct evolution under the imperfect controls. Formally, we may perform the 
expansion 



^2 rf2 



fMt); e] = fMty^ 0] + e^jMt); o] + |^/m[^W; o] + o{e 



(5) 



then by the condition that when e = the control must be error free, it is trivial to identify 
u^{t) = f^[u{t)]0]. Frequently, it is sufficient to consider models which are linear in the 
parameter e. In this case we introduce the shorthand notation 5u^{t) = ^/^[w(t);0] and 
the corresponding vector 5u{t) to represent the first-order deformation of the controls, so 
that the imperfect controls take the form v^(t) = Ufj,{t) + eSUfj,(t). Figure 2 illustrates two 
common error models: constant offsets from the ideal control value (/^[■u(t); e] = M^(t) + e), 
and errors in the control amplitude (/^[w(t); e] = (1 + e)'U^(t)). 

A natural question to ask is what effect unknown errors have on the evolution of the 
system. It is obvious that imperfect pulses make accurate manipulation of a quantum state 
difficult. One may be surprised to find that for some cases, the effects of errors on the 
controls may be systematically removed, without knowledge of the amplitude e. The method 
we describe involves implementing a compensating composite pulse sequence which is robust 
against distortion of the controls by a particular error model. As an example, consider a case 
where an experimentalist would like to approximate a the target unitary Ut = U{u{t); r, 0), 
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where at least one of the controls is influenced by a systematic error. The target operation 
may be simulated up to C(e") if there exists a set of control functions u{t) such that, 

V{u{ty, r, 0) = U{u{t) + e5u{ty, r, 0) = f/(w(t); r, 0) + 0(e"+^). (6) 

There are many (infinite in most cases) sets of control functions u{t) which implement a 
target unitary transformation, however only a small subset of possible control functions are 
robust to distortion by a particular systematic error. If robust controls can be found, then 
pulse sequence V{u(t); r, 0) can be applied in the place of Ut, and the leading order terms of 
the offset e6u(t) are suppressed. A sequence with these properties is called a compensating 
pulse sequence, and may be thought of as a set of control functions which are optimized to 
remove the effect of leading-order terms of unknown systematic errors in the controls. By 
construction, sequences of this form are fully compensating (Class A). 

When an experimentalist implements a compensating pulse sequence they attempt to 
apply the ideal operations U{u(t)]T,0) ignorant of the amplitude e of a systematic error. 
However, the operations are not ideal and to emphasize this we introduce V{u{t);T,0) to 
represent the imperfect propagators that is actually implemented when [/(■u(t); r, 0) is at- 
tempted. The functional dependence of the error in V{u{t);T,0) is not explicitly written 
allowing us to study different error models with the same pulse sequence. 

Let us consider the dynamics of the system under the interaction frame Hamiltonian 
H^it) = ^^e6Uf,{t)Hl^{t), where Hf^it) = W{u{t');t,0)H^U{u{t'y,t,0) are the control 
Hamiltonians in the interaction frame. In this picture, H\t) is regarded as a perturbation, 
and we associate the propagator U\Su{t); r, 0) as the particular solution to the interaction 
picture Schrodinger equation over the interval < t < r. Hence 

V{u{ty, r, 0) = U{u{ty T, 0)U\e6u{ty T, 0), (7) 

and from Eq. (6) 

f/^(e5w(t);r,0) = l + C(e"+^). (8) 

Quite generally, when a fully compensating pulse sequence is transformed into the inter- 
action frame the resulting propagator must approximate the identity operation [31]. The 
techniques for constructing compensating pulse sequences discussed in the present article 
rely on performing a series expansion by powers of e for the interaction frame propaga- 
tor U\e5u{t);T,0), then choosing a set of controls which remove the leading terms of the 
distortion, and finally transforming back out of the interaction frame. 

2.2 NMR spectroscopy as a model control system 

In this section, we will apply the ideas developed thus far to a model one-qubit NMR quantum 
computer [32, 33], which serves as a relevant example of a system where coherent control 
is possible. In section 5.1 multi-qubit operations are considered. There are many possible 
physical implementations for a qubit, oftentimes based on a two- level subsystem of a larger 
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Figure 3: Simplified diagram of an NMR spectrometer. Control over an ensemble of nuclear 
spins is applied by a radiofrequency (rf) field Bi{t), where the Rabi frequency fi{t), field 
detuning A(t), and phase 0(t) are control parameters. Unknown systematic control errors 
(e.g. poor intensity control, detuning errors) may be present. 



Hilbert space. In this case, the qubit is defined on the angular momentum states of a spin-1/2 
nucleus. 

Consider an ensemble of spin 1 = 1/2 nuclei undergoing Larmor precession under a static 
magnetic field Bo oriented along the z axis. The static field Bq induces a net magnetization 
among the nuclei. A transverse radiofrequency (rf) field near nuclear resonance Bi{t) = 
Bi cos{uJit — 0) is applied in the x-y plane. A simplified diagram of an NMR spectrometer 
is provided in figure 3. The analysis is simplified by assuming that individual nuclei in the 
ensemble are decoupled by the rapid tumbling of spins in the sample [25] . After transforming 
into the rotating frame, the Hamiltonian for a single spin may be written as 



H{t) = A{t)H, + n{t) ( cos(0(t))iy, + sin(0(t))i7y) , 



(9) 



where 



H. 



1 /o r 

2 u oj' 



H^ 



2 Vi 







H. 



1/1 
2 lo -1 



and the counter-rotating terms have been neglected under the rotating-wave approximation. 
Manipulation of the rf field provides a convenient set of controls to guide the evolution of 
the spins; it is assumed that the rf field detuning A(t), phase (p(t), and the Rabi frequency 
Q(t) are each independently controllable and may suffer from independent systematic er- 
rors. Let the controls u{t) denote the vector (fi(t) cos(0(t)), fi(t) sin(0(t)), A(t)). Then the 
Hamiltonian may be written as H{t) = u{t) ■ H, and the resulting Schrodinger equation for 
the propagator is of the form of Eq. (1). Thus, the coherent spin dynamics in magnetic 
resonance spectroscopy can be reformulated in terms of a problem in quantum control. 
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2.2.1 Error models 

In practice, systematic errors in the controls caused by instrumental limitations prohibit 
the application of perfect pulse propagators. We consider several models for errors in the 
controls of a one-qubit NMR quantum computer. Often compensating pulse sequences are 
well optimized for one type of error, but provide no advantage against a different error model. 

Amplitude errors: An amplitude error arises from slow systematic variation in the am- 
plitude of the rf field, resulting in a small offset in the applied Rabi frequency. Let Q 
represent the ideal Rabi frequency and Q' represent the offset. From the form of the control 
vector u{t) it follows that in the presence of the error, the Ux{t) and %(t) controls are dis- 
torted, i.e. they are replaced by their imperfect counterparts v^/yit) = u^/yit) + eA^u^/yit), 
where 6ux/y{t) = Ux/y{t) is proportional to the ideal control value and the error parameter 
eA = fi'/^ < 1 is the relative amplitude of the field offset. The imperfect pulses take the 
form V{u{t)) = U{u{t){l + eA)). 

Pulse length errors: A pulse length error is a systematic error in the duration of indi- 
vidual pulses, perhaps due to an offset in the reference oscillator frequency. The imperfect 
propagator takes the form V{u(t);to + At, to) = U{u(t);tQ + At + St, to), where At is the 
ideal pulse length and St is the unknown timing error. In some cases, errors on the clock 
may be rewritten in terms of equivalent errors on the control functions. For simplicity, we 
restrict ourselves to square pulses, where the controls u{t) are constant and the imperfect 
propagator may be rewritten as V{u; tf, ti) = U^u + erSu; t^ + At, to), where again Su = u is 
proportional to the ideal control and now e^ = St/ At, similar to the result for an amplitude 
error. Amplitude and pulse length errors in square pulses are similar in some senses, though 
they arise from distinct physical processes, since they both act as errors in the angle of the 
applied rotation. 

Addressing errors: Individual spins among the ensemble may experience slightly different 
Rabi couplings owing to the spatial variation in the strength of the control field. In some 
applications this variation is exploited to yield spatially localized coherent operations, such 
as in magnetic resonance imaging, and in addressing single atoms in optical lattices or single 
ions in ion trap experiments [34, 17]. In fact, many proposed scalable architectures for 
quantum processors rely on this effect to discriminate between qubits. In these cases, it 
is important to distinguish between the evolution applied to the addressed spins and the 
evolution of spins outside of the addressed region, where ideally no operation is applied. 

Consider an experiment where the ensemble of spins is divided into a high field, Q, and 
low field, Q' < Q, region. We assume that the control of the addressed (high-field) spins 
is perfect (i.e., for addressed spins V{u(t)) = U{u(t))); however, the spins in the low-field 
region experience an undesired correlated rotation V{u{t)) = U{eNu{t)), where cat = Q! /Vt 
(the subscript A^ denotes neighboring spins). In many cases, the imperfect pulses on the 
unaddressed qubits may be regarded as very small rotations. From a mathematical point of 
view, sequences composed of these rotations are more easily attacked using the expansion 
techniques that will be developed in section 3.2. 

Detuning errors: Systematic errors may also arise in the control of the frequency of 
the rf field. Consider the case where an experimentalist attempts to perform an operation 
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at a particular field tuning A(t), however a slow unknown frequency drift 6{t) is present. 
From the form of the NMR controls, it follows that the error distorts the H^ component 
to v^it) = u^{t) + e£)6uz, where u^ = A(t), and eodu^it) = 6{t). Unlike the other models 
discussed thus far, the detuning error applies a shift in the effective rotation axis, rotating it 
in the direction of the H^ axis. Therefore, the ideal propagator and the imperfect counterpart 
do not commute in general. 

2.3 Binary operations on unitary operators 

In our discussion of the properties of compensating pulse sequences, it will be necessary to 
study the structure of the control Hamiltonians {-ff/x}, and to calculate the accuracy of gates. 
Here we discuss the Hilbert-Schmidt product and the fidelity measure. 

Hubert- Schmidt inner product: The Hilbert-Schmidt inner product (also known as the 
Frobenius product) is a natural extension of the vector inner product over the field of complex 
numbers to matrices with complex coefficients. Let U and V he n x n matrices with entries 
Uij and Vij respectively. The Hilbert-Schmidt inner product {U, V) is defined as, 

n n 

{u,v) = J2T.u;^v^^ = ^'(u^^y (10) 

This inner product has several properties that closely mirror the inner product for complex 
vectors, namely {U,V) = {V,U)*, and \{U,V)\ < \\U\\us H'^^Hhs, where the Hilbert-Schmidt 
norm ||6'"||hs = ^/{U, U) also satisfies a corresponding triangle inequality ||f/ + ^||hs < 
II^IIhs + II^IIhs- This strong correspondence to Euclidean vector spaces will be exploited 
in section 3.1, where Lie-algebraic techniques are employed to construct compensating pulse 
sequences. 

Fidelity: A useful measure for evaluating the effects of a systematic control error is the 
operational fidelity, defined as 

J^{U, V) = min ^y{tlj\WV\tlj){ilj\Vm\ilj). (11) 

For any two unitary matrices in the group SU(2), the fidelity may be written as J^{U, V) = 
\{U,V)\/2. In this paper, we show error improvement graphically by plotting the infidelity, 
1-T{U,V). 

3 Group theoretic techniques for sequence design 

In the present work, we briefly discuss several aspects of the theory of Lie groups which are 
useful for constructing compensating pulse sequences, with emphasis on conceptual clarity 
over mathematical rigor. The interested reader is referred to Refs. [10] and [35] for additional 
details and a rigorous treatment of the subject. Quantum control theory and the related 
study of continuous transformation groups is a rich and active subject, with applications in 
chemistry [36, 37, 38], and other applications requiring the precise manipulation of quantum 
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states [39, 40]. We now turn our attention to continuous groups of quantum transformations, 
specifically transformations which induce qubit rotations in quantum computing. 

3.1 Lie groups and algebras 

Consider the family of unitary operators that are solutions to the control equation 

U{t) = Y,Mt)H,U{t), f/(0) = l, (12) 

where H^^ = —iH^ are skew-symmetrized Hamiltonians with the added condition that the 
(possibly infinite) set of control Hamiltonians {H} is closed under the commutation operation 
(in the sense that for all if^, Hy G {i?}, a[H^,Hy\ G {H}, where a G R ). This extra 
condition does not exclude any of the "physical" operators which are generated by a subset 
of {-ff}, since these operators correspond to solutions where the remaining control functions 
u^{t) have been set to zero. Observe that these solutions form a representation of a group, 
here denoted as G, since the following properties are satisfied: for all solutions f/i and f/2, the 
product U1U2 is also a solution (see section 3.2); the associative property is preserved (i.e. 
Ui{U2U^) = {UiU2)Us); the identity 1 is a valid solution; and for all solutions Ui, the inverse 
Ul is also a valid solution. Moreover, the group forms a continuous differentiable manifold, 
parametrized by the control functions. A continuous group which is also a differentiable 
manifold with analytic group multiplication and group inverse operations is called a Lie 
group. We identify the group G of solutions as a Lie group, and note that differentiation of 
the elements U E G is well defined by nature of Eq. (12). 

We now turn our attention to elements of G in the neighborhood of the identity element, 
that is, infinitesimal unitary operations. In analogy with differentiation on Euclidean spaces, 
note that for any U{t) G G one may find a family of tangent curves at f/(0) = 1, 



dU{t) 



dt 



t=o 



j2Mm,- (13) 



The set of skew-symmetrized Hamiltonians {H} and the field of real numbers R (correspond- 
ing to the allowed values for the components ^^^(O)) form a linear vector space under matrix 
addition and the Hilbert-Schmidt inner product [35]. Here, the Hamiltonians H take the 
place of Euclidean vectors and span the tangent space of G at the identity, denoted by TjG. 
A homomorphism exists between the control functions u{0) and vectors in TjG. On this 
space is defined the binary Lie bracket operation between two vectors [Hfj,,Hy], which for 
our purposes is synonymous with the operator commutator. A vector space which is is closed 
under the Lie bracket is an example of a Lie algebra. By construction, the set {H} is closed 
under commutation and therefore forms a Lie algebra, here denoted as q, corresponding to 
the Lie group G. We note that in general the Hamiltonians may be linearly dependent; 
however, an orthogonal basis under the Hilbert-Schmidt product may be generated using an 
orthogonalization algorithm. 
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The power of most Lie algebraic techniques rehes on the mapping between group elements 
which act on a manifold (such as the manifold of rotations on a Bloch sphere) to elements 
in a Lie algebra which are members of a vector space. In the groups we study here, the 
mapping is provided by the exponential function G = e^ (i.e., every element [/ G G may be 
written as [/ = e^ where g & q). The object of this method is to study the properties of 
composite pulse sequences, which are products of members of a Lie group, in terms of vector 
operations on the associated Lie algebra. 

3.1.1 The spinor rotation group SU(2) 

As a relevant example, consider the group of single-qubit operations generated by the NMR 
Hamiltonian Eq. (9). This is a representation of the special unitary group SU(2). The 
skew-symmetrized control Hamiltonians are closed under the Lie bracket and thus form a 
representation of the Lie algebra 5u(2) = span{— ii^x? ~^Hy, —iH^} ^. Therefore, any element 
in ?7 G SU(2) may be written as f/ = e"'*'"'-^, where —ituH G su(2) may now be interpreted 
as a vector on the the Lie algebra. Furthermore, since (— iif^, —iH^) = 5^^y/2, the spin 
operators form an orthogonal basis for the algebra. Topologically SU(2) is compact and is 
homomorphic to rotations of the 2-sphere, the group SO (3) [35]; there are exactly two group 
elements U and —U G SU(2) which map to the same rotation. 

3.2 Baker-Campbell-Hausdorff and Magnus formulas 

A composite pulse sequence may be studied from the perspective of successive products be- 
tween elements of a Lie group. In the following analysis, it will be useful to relate the product 
of two members of a Lie group to vector operations on the Lie algebra. The relationship 
allows us to map pulse sequences to effective Hamiltonians. This correspondence is provided 
by the Baker-Campbell-Hausdorff (BCH) formula [35, 41], which relates group products to a 
series expansion in the Lie algebra. For rapid convergence, it is most convenient to consider 
products of infinitesimal unitary operations, that is, operations of the form e*^^, where f? G 
and e < 1 is a real expansion parameter. We assume e is sufficiently small to guarantee 
that the group product of propagators always lies within the radius of convergence for the 
expansion [20]. Let Ui = e'^^^ and f/2 = 6*^^^ be members of a Lie group G, where the 
Hamiltonians Hi, H2 ^ Q are members of the associated algebra. The BCH representation 
for the product U1U2 = U3 involves the calculation of an effective Hamiltonian i^s G by 
the expansion 

U, = exp(^3) = exp ( f^ e"F„ ) , (14) 



^The operation span denotes all linear combinations with real coefficients 



3 GROUP THEORETIC TECHNIQUES FOR SEQUENCE DESIGN 13 

where the terms 



eFi = e{Hi + H. 



2 



e2 



e'F2 = -[H,,H2] 



e^Fs = ^([i^i,[^i,^2]] + [i^2,[i^2,i^l] 



are calculated from Hi, H2, and nested commutators of elements of the Lie algebra. A com- 
binatoric formula found by Dykin [42] exists to calculate Fn for arbitrary n. The expansion 
may be truncated once a desired level of accuracy is reached. In principle, group products of 
arbitrarily length may be approximated to arbitrary accuracy using BCH formulas; however, 
these formulas rapidly become unwieldy and difficult to use without the aid of a computer. 
The BCH expansion is most useful for sequences of square pulses, where for each pulse the 
Hamiltonian is time independent. 

The Magnus expansion: A related expansion developed by Magnus [19] may be used to 
compute the propagator generated by a general time-dependent Hamiltonian. The solution 
to a control equation (e.g. U{t) = eH{t)U(t), where eH(t) = ^ e5u^{t)H^) over the interval 
ti < t < tf may be written as the power series 

U{eSu; tf, U) = exp I ^ e"n„(t/, U) J , (15) 

where the first few expansion terms are, 

efii(t/,ti) = e / dtH{t) 
Jti 

c2 /•</ /•< 



2 ft f PL 

e^n^itf^U) = '- dt dt'[H{t),H{t')] 

e'ns{tf,U) = ^ l^'dt f dt' f dt"{[H{t),[H{t'),H{t")]] + [H{t"),[H{t'),Hm)- 



Again as a matter of notational convenience, we drop the time interval labels {tf,ti) on 
the expansion terms when there is no risk of confusion. Formulas for higher order terms 
may be found in [43]. The BCH and Magnus expansions are in fact intimately related; when 
considering piecewise-constant controls the techniques are equivalent. We refer the interested 
reader to [20] for further details regarding both the Magnus expansion and BCH formulas. 
The BCH and Magnus expansions are very well known in composite pulse literature, and 
techniques that utilize these expansions are collectively referred to as average Hamiltonian 
theory. A variant of this technique, pioneered by Waugh [44, 45], has been a mainstay of 
composite pulse design in the NMR community for decades. In some formalisms, the BCH 
expansion of the product of two propagators is interpreted as a power series in the rotation 
axis and angle. For the study of composite single-qubit rotations, this picture is extremely 
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Figure 4: Vector paths on a Lie algebra q may be used to represent a pulse sequence, a) 
The BCH expansion relates group multiplication of several square pulses to vector addition 
on 0. b) A shaped pulse is represented by a vector curve on q, parametrically defined by the 
control functions. Both pulse sequences form a closed loop on q and therefore the first-order 
terms equal zero. 

useful as it allows rotations on the sphere to guide the mathematics. However, this picture 
of composite rotations can not be generalized to more complex groups, such as the group of 
n-qubit operations SU(2"). 

In this work, we emphasize a Lie algebraic interpretation of these methods, which also 
leads to a second geometric picture for the terms of the expansions. Observe that the first- 
order terms Fi and fli may be regarded as simple vector sums on the Lie algebra, i.e., the sum 
of Hi and H2 in the BCH expansion, and the sum of each of the infinitesimal vectors H{t)dt 
in the Magnus expansion. In an analogous way, one may interpret the higher order terms 
as the addition of successively smaller vectors on q. From this insight, one may construct 
composite sequences which simulate a target unitary from geometric considerations on the 
Lie algebra. 

3.2.1 A method for studying compensation sequences 

At this point, it is useful to introduce the general method that will be used to study compen- 
sating pulse sequences. Recall from section 2.1 that in the presence of an unknown systematic 
control error, the ideal control functions are deformed into imperfect analogues. The imper- 
fect propagator may be decomposed as V{u{t)]T,Q) = U{u(t)]T,0)U\eSu{t);T,0), where 
U\e5u{t); r, 0) represents the portion of the evolution produced by the systematic distortion 
of the ideal controls. Provided that the displacements e5u{t) are sufficiently small relative to 
the ideal controls, we may perform a Magnus expansion for the interaction frame propagator 



U\t8u{t)] T, 0) = exp I ^ e"fin(r, 0) 



(16) 



vn=l 
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where the integrations in the expansion terms are performed in the appropriate frame. For 
example, the first-order term is 

efii(r, 0) = -ie f dt Su{t) ■ H\t), (17) 

Jo 

where the components Hl^it) = U'^{u(t');t,0)H^U{u(t');t,0) are the interaction frame con- 
trol Hamiltonians. The reader should recall from Eq. (8) that if the controls u{t) form an nth 
order compensating pulse sequence, then the interaction frame propagator U\e6u{t);T, 0) = 
1 + 0{e^~^^) must approximate the identity to sufficient accuracy. It immediately follows 
that this condition is satisfied for any e when the leading n-many Magnus expansion terms 
Qi{t, 0), r22(r, 0), . . . , r2„(r, 0) over the pulse interval < t < r simultaneously equal zero. 

This condition may also be understood in terms of geometric properties of vector paths 
on the Lie algebra. For example, consider the first-order Magnus expansion term for the 
interaction frame propagator and observe that —ie5u{t) ■ H\t) may be regarded as a vector 
path on the dynamical Lie algebra. From the condition Qi{t,0) = and Eq. (17), we see 
that the path must form a closed cycle in order for the first-order term to be eliminated (see 
figure 4). The elimination of higher-order expansion terms will place additional geometric 
constraints on the path which will depend on the structure of the Lie algebra (i.e. the 
commutators between paths on the algebra). In the case of piecewise constant control 
functions, the resulting propagator may be understood as a sequence of square pulses. It 
is clear that in this case the Magnus series reduces to a BCH expansion for the total pulse 
propagator, and on the Lie algebra the corresponding path forms a closed polygon. 

This Lie theoretic method is a useful tool in determining whether a control function u{t) 
is also a compensating sequence; we may directly calculate the interaction frame Magnus 
expansion terms in a given error model and show that they equal zero. The inverse problem 
(i.e. solving for control functions) is typically much more difficult. In general, the interaction 
Hamiltonians HUt) are highly nonlinear functions of the ideal controls it(t), which impedes 
several analytical solution methods. However, we note that in several special cases the 
problem is considerably simplified, such as when the ideal operation is to perform the identity 
and the interaction and Schrodinger pictures are equivalent. 

3.3 Decompositions and approximation methods 

Several useful techniques in sequence design involve decompositions that may be understood 
in terms of the structure of a Lie group and its corresponding algebra. In this section we 
discuss several important methods and identities. 

3.3.1 Basic building operations 

Given a limited set of controls {Hi,H2} that generate the algebra g, one may produce any 
unitary operation in the corresponding Lie group G = e^ using only two identities. The first 
identity, the Lie- Trotter formula [46], describes how to produce a unitary generated by the 
sum of two non-commuting control operators. Using the BCH formula one may compute that 
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QHi/n^H2/n _ Q(Hi+H2)/n _|_ Q(^^]j-^^ jj^y^'^y In terms of physical pulses, this corresponds to 
dividing the propagators e^^ and e^^ into n equal intervals to produce the propagators e^^/" 
and e^^/". Suppose we perform n such successive products, so that the resulting propagator 
is 

LH^/n^H2/ny ^ ^H.+H^ ^ 0{[Hi, H2]/ u). (18) 

Although the Haniiltonians Hi and H2 do not commute in general, we may approximate 
U = e^^~^^^ to arbitrary accuracy by dividing the evolution into ra-many time intervals and 
using the construction Eq. (18). By extension, it follows that any unitary generated by a 
Hamiltonian in the Lie algebra subspace spa.n{Hi, H2} may be approximated to arbitrary 
accuracy using a Trotter sequence. A number of improved sequences were developed by 
Suzuki [47] that remove errors to higher commutators and scale more strongly with n. The 
Trotter-Suzuki formulas may be used to eliminate successively higher order errors at the cost 
of increased operation time [48] . 

The second identity, which we refer to as the balanced group commutator, enables the 
synthesis of a unitary generated by the Lie bracket [Hi,H2]. Again the BCH formula may 
be used to show e^i/"e^2/ng-i?i/ng-H2/n ^ ^iHum/n^ + Qi^^Hi + H2, [Hi,H2]]/n^). If we 
now consider n^-many successive balanced group commutator constructions, the resulting 
propagator is 

/ - ~ - ~ \ Tl ~ _ ^ ~ ~ ~ 

(^e^l/"e^2/ng-Hi/ng-H2/nj ^ ^[H^H^] ^ (>^Q^^ ^ ^^^ [^^^ H2]]/n). (19) 

Then, as in the case of the Trotter formula, we may approximate U = e^^^'^^^ to arbitrary 
accuracy by increasing the number of intervals n. Since by assumption the entire Lie algebra 
may be generated by nested Lie brackets between the Hamiltonians Hi and if 2, this implies 
that any f/ G G may be produced by a combination of balanced group commutator and 
Trotter formulas. However, we emphasize that in almost all cases much more efficient con- 
structions exist. The balanced group commutator construction also forms the basis of the 
Solovay-Kitaev theorem [49] , an important result regarding the universality of a finite gate 
set in quantum computation. Later, we will study the Solovay-Kitaev method (see 4.1.1) 
which produces compensating sequences of arbitrary accuracy by using a balanced group 
commutator. Specifically we will use the formula 

exp(^ie'=) exp(^2e') exp(-^ie'=) exp(-^2e') = exp([^i, ^2]e'+') + C(e*^+'+^), (20) 

where the parameter e < 1 will represent the strength of a systematic error. 

3.3.2 Euler decomposition 

In section 3.1.1 it was shown that any one-qubit operation U G SU(2) may be written in the 
form U = exp{—itu ■ H). It is well known that an alternative representation exists, namely 
the Euler decomposition 

U = exp(— iaaifx) exp(— ia2-ffy) exp(— iaiifx), (21) 
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which is given by sequential rotations by the angles {01,02,0:3} about the H^, Hy, and 
Hyr axes of the Bloch sphere. The Euler decomposition gives the form of a pulse sequence 
that produces any arbitrary one-qubit gate. The Euler decomposition also gives a method 
of producing rotations generated by a Hamiltonian outside of direct control. For example, 
at perfect resonance U = exp{—i6Hz) cannot be directly produced; however, the Euler 
decomposition U = exp(— i|iJx) exp(— i6'iJy) exp(i|_fZ'x) may be implemented. 

3.3.3 Cartan decomposition 

Let be a semi-simple Lie algebra that may be decomposed into two subspaces g = { © m, 
m = {"*" satisfying the commutation relations, 

% I] C I, [m, I] C m, [m, m] C t (22) 

Such a decomposition is called a Cartan decomposition of q [10]. Suppose for the moment 
there exists a subalgebra a of g which is in a subspace of m. Since a is an algebra of its 
own right, it is closed under the Lie bracket [a, a] C a. However, note C m implies that 
[CI7 1] ^ [tTi, m] C t. Since the subspaces t and m are mutually orthogonal, then [a, a] = {0} 
and the subalgebra a must be abelian. A maximal abelian subalgebra a C m for a Cartan 
decomposition pair {t,m) is called a Cartan subalgebra [10]. 

For brevity, we state without proof an important theorem regarding the decomposition 
of an operator in a group G with a Lie algebra admitting a Cartan decomposition. Consider 
a Lie algebra q with a Cartan subalgebra a corresponding to the decomposition pair (J^, m). 
Every U in the group G = e^ may be written in the form, 

U = K^AK^, (23) 

where Ki,K2 G e' and A G e". This is called the KAK Cartan decomposition for the 
group G. The interested reader is referred to [35] for additional details regarding the KAK 
decomposition. 

As a relevant example, here we show how the Euler decomposition for a propagator U G 
SU(2) is a special case of a KAK decomposition. The algebra is spanned by the orthogonal 
basis matrices su(2) = span{— iiJx, ~^Hy, —iH^}. Observe that t = span{— iifx} and m = 
span{— iify, —iH^} form a Cartan decomposition for 5u(2). The maximal abelian subalgebra 
of m is one-dimensional; we choose a = span{— iify} although the choice a' = span{— iif^} 
would serve just as well (i.e. the different axis conventions of the Euler decomposition 
differ in the choice of a maximal abelian subalgebra). Then by Eq. (23), every element 
U G SU(2) may be expressed in the form U = exp(— ioaifx) exp(— io2-ffy) exp(— ioiilfx), 
where the parameters Oj are real. This is a restatement of Eq. (21), thus completing the 
proof. The KAK decomposition is an existence theorem, and does not provide a direct 
method for the calculation of the required rotation angles Oj. 

The Cartan decomposition has important implications for universality. For instance, if 
one may generate any unitary operation over the subgroups e' and e", then any gate in the 
larger group e^ may be produced. Similarly, if compensation sequences exist for operations 
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in these subgroups, then they may be combined in the KAK form to yield a compensating 
pulse sequence for an operation in the larger group. Another important application is the 
decomposition of large Lie groups into products of more simple ones. Of special interest to 
quantum computing is the inductive decomposition of n-qubit SU(2") gates into products of 
one-qubit SU(2) and two-qubit SU(4) rotations [50]. 

4 Composite pulse sequences on SU(2) 

In this section, we study pulse sequences which compensate single-qubit operations, which 
form a representation of the group SU(2). Our approach is to use Lie theoretic methods to 
study the effects of composite rotations. Topologically, sequences of infinitesimal rotations 
may be interpreted as paths in the neighborhood of the group identity. It is frequently 
easier to construct sequences of infinitesimal paths on the Lie algebra, and then map these 
sequences to the manifold of group operations. Several techniques will be used to construct 
compensating pulse sequences of arbitrary accuracy, including techniques that use Solovay- 
Kitaev methods and Trotter formulas. 

4.1 Solovay-Kitaev sequences 

The Solovay-Kitaev sequences, so named because the construction of the higher order se- 
quences involves the identity Eq. (20) used by Kitaev in his proof of universal control from 
finite gate sets [49] , are among the simplest families of fully-compensating composite pulse 
sequences. The SK sequences were first introduced by Brown, Harrow, and Chuang [51], and 
are designed to compensate pulse length, amplitude, and addressing errors using resonant 
pulses. Here we show that the SK family of sequences may be derived using the method 
outlined in section 3.2.1 

Narrowband behavior: Narrowband composite pulse sequences apply a spin rotation over 
only a narrow range of strengths of the control field [52]. Therefore, they are most suited 
for correcting addressing errors [34], i.e., situations where the spatial variation in the field 
strength is used to discriminate between spins in an ensemble. The operations on the ad- 
dressed qubits are assumed to be error-free, whereas on the unaddressed qubits the imperfect 
pulses take the form V{u(t)) = U(eNu{t)), where ejv < 1 is the systematic addressing er- 
ror amplitude (see section 2.2.1). Narrowband sequences are a means of applying a target 
operation Ut on the addressed spins while removing the leading effects of the operation on 
the unaddressed spins. For a composite pulse sequence to exhibit nth-order narrowband 
behavior, we require that two conditions must be satisfied: first that for the unaddressed 
qubits the sequence V{u{t)) approximates the identity up to 0{e'^^), and second, that on 
the targeted spins the desired operation is applied without error. 

From a mathematical perspective, systematic addressing errors are among the easiest to 
consider. Note that for the unaddressed qubits the ideal values for the controls u{t) is zero 
(ideally no operation takes place). This implies -ff^(t) = H^, and that the Schrodinger and 
interaction frames in Eq. (8) are identical. We may develop a sequence which corrects the 



4 COMPOSITE PULSE SEQ UENCES OiV S U (2) 19 

control distortion without the added comphcation of the passage into an interaction picture. 
For this reason, we first study narrowband sequences before considering other error models. 
We now show that a sequence with first-order narrowband properties may be constructed 
using the method described in section 3.2.1. We consider sequences composed of three square 
pulses (i.e., piecewise constant control functions) which produce the ideal propagation, 

3 

U{u{t)) = II Uk, Uk = exp{-itkUk ■ H). (24) 

fc=i 

On an addressed qubit, U{u{t)) is implemented, whereas on the unaddressed qubits V{u{t)) = 
U{eNu(t)) is applied. We may use either a BCH or Magnus expansion to compute the applied 
operation on the unaddressed qubit. From Eq. (17) and the error model fx/y(^) = ^NU^/yif)-, 
the first-order term is 

eAT^i = -itN^hu^i + t2U2 + tiUi) ■ H. (25) 

The spin operators —iH^ G {— ii/x, ~^Hy, —iHz} form an orthogonal basis for su(2), and the 
terms —iej^tkUk ■ H may be regarded as vectors on the dynamical Lie algebra. In order to 
eliminate the first-order Magnus term eArfii, the sum of the components must equal zero, 
that is, the vectors must form a closed triangular path. Such paths may be found using 
elementary geometric methods. 

It is important at this point to allow experimental considerations to place constraints on 
the sequences under study. For instance, in many cases it is desirable to perform coherent 
operations at resonance, a condition which forces the vectors to lie in the Hy^-Hy plane. One 
possible choice is 

— leNtiUi ■ H = —iejqOH^ 

-ieNt2U2-H = -ieAr27r(cos0sKi-f^x + sin0sKi-f^y) 

-ieNtsUs-H = -ieAr27r(cos0sKi-f^x-sin0sKi-f^y), (26) 

where the phase 0ski = arccos(— 6'/47r) is selected so that e^- Ylik^kUk = 0, and therefore the 
first-order expansion term e^^i = is eliminated. Figure 5a is a diagram of these vectors on 
5u(2), where the sequence may be represented as a closed isosceles triangle with one segment 
aligned on the —iH^ axis. 

For clarity, we use the notation R{d, (p) = exp{— i6 [cos (f)Hy + sin (pHy)) to represent prop- 
agators which induce rotations about an axis in the H^-Hy plane,. Observe that for resonant 
square-pulse operations in SU(2), R{6,(j)) = U{uk;tk,0), where Uk = f2(cos0, sin0, 0) and 
9 = tk\uk\. We also define the corresponding imperfect propagator M{9, 0) = V{u; tk, 0), and 
recall that the imperfect propagators on the unaddressed qubits M{9, 0) = R{9eN, 0) and 
addressed qubits M{9, 0) = R{9, 0) have different implied dependencies on the systematic 
error. Then combining Eq. (24) and Eq. (26), the propagator may be written as 

U{eNu{t)) = R{27ieN, -^sKi)R{27reN, (psKi)R{9eN, 0) = 1 + 0{e%) 
MsKi(^,0) = M(27r,-0sKi)M(27r,0sKi)M(^,O), (27) 
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Figure 5: a) Vector path followed by SKI on the Lie algebra, b) Trajectory of an unaddressed 
spin during an SKI sequence, using imperfect rotations of the form M{9,(l)) = R{9eiy,(f)). 
c) SKI correcting an amplitude error, using imperfect M{9,(f)) = R{9{1 + e^),0). In these 
plots eN = €a = 0.2. 



for the unaddressed qubits, and as 

U{u{t)) = R{2tt,^ 
MsKi{9,0) = M(27r, 



bsKl)Ri27T,(l)sKl)Ri9,0) 

0SKi)M(27r,0sKi)M(^,O), 



(28) 



for addressed qubits. This is the first-order Solovay-Kitaev sequence, here denoted as SKI. 
On the addressed qubits the rotations i?(27r,0sKi) = -R(27r, — 0ski) = —1 are resolutions 
of the identity and the effect of the sequence is to apply the target unitary Ut = R{9,0), 
whereas on the unaddressed spins the sequence applies 1 + 0{e%) and the leading first-order 
rotation of the unaddressed spins is eliminated. Therefore, SKI satisfies the conditions for 
narrowband behavior. When the sequence SKI is used in the place of the simple rotation 
R{9,0), the discrimination between addressed and unaddressed spins is enhanced. 

Broadband behavior: Broadband composite pulses apply a spin rotation over a large range 
of strengths of the control field, and are best suited for correcting systematic amplitude 
and pulse-length errors that correspond to systematic over/under-rotations during qubit 
manipulations. Broadband sequences are a means of applying a target operation in the 
presence of inaccurate field strengths or pulse durations. We require for a sequence to 
exhibit nth-order broadband behavior, the effect of the sequence is to approximate the 
target operation Ut up to 0{e'\) in the case of amplitude errors, or up to 0{e^) in the case 
of pulse-length errors. As before, it is convenient to consider sequences comprised of resonant 
square pulses. In this case the amplitude and pulse length error models are in some sense 
equivalent since they both apply a proportional distortion to the u^it) and Uy(t) controls. 
In the following discussion we will explicitly use the amplitude error model, where imperfect 
rotations take the form M{9, 0) = R{9{1 + e^), 0). 

We now show how a Solovay-Kitaev sequence with first-order broadband properties may 
be derived. Although the method presented in section 3.2.1 may be applied, here a more 
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direct technique is used. Our strategy is to construct a pulse sequence entirely out of imper- 
fect rotations M{e,(j)) = R{6{1 + e^),*/*) = R{9eA,(j))R{0A) using Eq. (27) as a template. 
Explicitly, this is achieved with the matrix manipulations 

R{9,0) + Oiel) = R{27reA,-(psKi)Ri27reA,(pSKi)Ri0eA,0)Rie,0) 

= Ri2TTeA, -0sKi)i?(27r, -0sKi)i?(27reA, 0sKi)i?(27r, 0sKi)i?(^e^, 0)i?(^, 0), 

where we have right-multiplied Eq. (27) by R{9, 0) (after substituting eA for e^v) and inserted 
the identity R{27i , (psKi) = -R(2vr, — 0ski) = —1- Then, by combining rotations about the 
same axis, one obtains the result, 

Uiil + eA)u{t)) = /?(27r(l + e^),-0sKi)i?(27r(l + eA),0SKi)i?(^(l + eA),O) 

MsKi(^,0) = M(27r,-0sKi)M(27r,0sKi)M(e,O)=i?(^,O) + O(ei) (29) 

where again 0ski = arccos(— ^/47r). In the presence of unknown amplitude errors, the 
sequence reduces the effect of the error while applying the effective rotation R{6, 0). Observe 
that in terms of imperfect rotations this is the same SKI sequence derived earlier. SKI has 
both narrowband and broadband behavior, and the sequence may correct both addressing 
and amplitude errors simultaneously. Such a sequence is called a passband sequence. 

Generalization to arbitrary gates in SU(2).- The sequence SKI is designed to compen- 
sate single qubit rotations about the H^ axis. If this sequence is to be useful in quantum 
computation, it must be generalized so that any single qubit rotation may be corrected. 
One method involves transforming the sequence by a similarity transformation of the pulse 
propagators. As an example, suppose we require an SKI sequence that performs the rotation 
R{9, 0) = exp{—i(j)Hz)R{9, 0) exp(i0i7z) on the addressed spins. Similarity transformation of 
pulses under exp(— i^if^) represent a phase advance in the rotating frame. It is evident that 
the transformed sequence 

MsKiiO, 0) = M(27r, - 0sKi)M(27r, + 0ski)M(^, 0) (30) 

performs the desired compensated rotation. In this manner, a compensating pulse sequence 
for any target operation Ut € SU(2) may be constructed. One may solve for the operation T 
that performs the planar rotation Tt/j'T"'' = R{9,0), where 9 = || log[/T||Hs/||-f^x||HS- Then 
the transformed sequence T'''Mski(6', 0)T performs a first-order compensated Ut operation. 
Similarity transformations of pulse sequence propagators can be extremely useful, and are 
frequently applied in composite pulse sequences. 

Alternative methods exist for generating an arbitrary compensated rotations. Recall from 
section 3.3.2 that any operation Ut € SU(2) may be expressed in terms of a Euler decompo- 
sition Ut = -R(«3, 0)R{a2, 'ir/2)R{ai, 0). An experimentalist may apply Ut by implementing 
each of the Euler rotations in sequence, however in the presence of an unknown systematic 
error each of the applied rotations is imperfect and the fidelity of the applied gate is reduced. 
The error may be compensated by replacing each imperfect pulse with a compensating pulse 
sequence. For example, the sequence 

MsKi(a3, 0)MsKi(a2, 7r/2)MsKi(ai, 0) = Ut + 0{e\) (31) 
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compensates amplitude errors to first-order by implementing SKI sequences for each of the 
rotations in an Euler decomposition for Ut- This construction is an example of pulse sequence 
concatenation. By concatenating two independent pulse sequences, it is sometimes possible 
to produce a sequence with properties inherited from each parent sequence. 

4.1.1 Arbitrarily accurate SK sequences 

In this section, we discuss the Solovay-Kitaev method for constructing arbitrarily accurate 
composite pulse sequences, which may be used to systematically improve the performance of 
an initial seed sequence. The Lie algebraic picture is particularly helpful in the description 
of the algorithm. The method is quite general, and can be used on sequences other than 
SKI. 

Suppose that we have an nth-order compensating pulse sequence, here denoted as Wn- 
The problem we consider is the identification of a unitary operator An+i such that Wn+i = 
An+iWn = Ut + C(e""'"^), where Wn+i is an (n + l)th-order sequence. Assume for now 
that such an operator exists, and consider that in the presence of systematic errors, the 
application of the correction gate An+i is imperfect. However, if it is possible to implement 
a compensating pulse sequence Bn+i, which is an C(e""'"^) approximation of An+i, then it is 
still possible to construct an (n + l)th-order sequence, . 

Wn+i = B,,+,Wn = Ut + 0{e''+^). (32) 

We may then continue constructing pulse sequences of increasing accuracy in this fashion if 
there exists a family of operators {An+i, An+2, ^n+s, ■ ■ ■ Am} and a corresponding family of 
pulse sequences {-B„+i, -Bn+25 -Bn+s, ■ ■ ■ -Bm} which implement the operators to the required 
accuracy. This immediately suggests an inductive construction for the sequence Wm, 

Wm = BmBm-l- ■ ■ Bn+-iBn+2Bn+lWn = UT + 0{e™' ). (33) 

This is the basis of the Solovay-Kitaev method [51]. To apply the method, we must first 
have a means of calculating the correction A„_|_i, and second we must find a compensating 
sequence -Bn+i robust to the systematic error model considered. 

We turn our attention to the calculation of the correction terms An+i- It is convenient to 
decompose the pulse sequence propagator using an interaction frame as Wn = UTU\e6u{t)), 
where Ut = U{u{t)) is the target gate. In the spirit of Eq. (16), a Magnus expansion for 
U^{e5u{t)) may be used. Observe that Wn = f/Texp(e"+in„+i) + C(e"+2). Then letting 
An+i = Ut exp{—e'^~^^iln+i)UT it may be verified using the BCH formula that An+iWn = 
Ut + 0{e"'~^'^). This result may also be interpreted in terms of vector displacements on the 
Lie algebra. In the interaction frame, e"'~^^Qn+i may be interpreted as a vector in su(2). 
Similarly, the infinitesimal rotation An+i corresponds to the vector — e"^^fi„+i, equal in 
magnitude and opposite in orientation. The first-order term of the BCH series corresponds 
to vector addition on the Lie algebra, and the 0{e"'~^^) terms cancel. 

What remains is to develop a compensating pulse sequence that implements An+i to 
the required accuracy under a given error model. Let us define Pjz{o:) = exp{—iae^ H^) + 
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0{e^'^^) and also the rotated analogues Pjxia) = exp(— iae-'ifx) + 0{e^~^^) and Pjy{a) = 
exp(— iae-'ify) + 0{e^^^). Frequently if one such Pj may be produced, then often the re- 
maining two may be produced by similarity transformation of the propagators or by using 
an Euler decomposition. At this point we use Eq. (20), used in the proof for the Solovay- 
Kitaev theorem, to construct relation 

PU-a)Piy{~^)PkA(^)Pey{/3) = Pjz{a/3), k + £ = j. (34) 

Continuing in this manner, each of the Pj's may be recursively decomposed into a product 
of first-order propagators Pik{(y) = exp(— iaei^k) + C'(e^)- Our strategy is to use Eq. (34) 
to implement P(^n+i)z{0 where ^ = ||f^n+i||Hs/||-f^z||HS- On the Lie algebra, this operation is 
represented by vector of length e""''^||i7„+i||Hs oriented along the —iH^ axis. Let T be the 
rotation which performs Ti7„+iT"'' = i^Hz (i.e., the operator which rotates Qn+i onto the 
iHz axis). Then, by similarity transformation under UtT"^ 

?7TTtP(„+i),(0T4 = f/Texp(-e"+if]„+i)4 + 0(e"+2) = b„.+i, (35) 

the sequence P(^n+i)z{0 '^^Y be transformed into precisely the required correction sequence 
needed for the Solovay-Kitaev method. There are two approaches for applying the transfor- 
mation by f/^T^: we may either calculate the transformed analogues of each of the pulses 
in P(n+i)z{0 ^^d then apply the transformed pulses, or we may, when possible, directly 
include the transformation pulses (or an estimate for them) as physically applied pulses in 
the sequence. The second approach is only viable when it is possible to generate accurate 
inverse operations TUj. [53]. 

Our construction will be complete once we have a method for generating the simple "pure- 
error" propagator Pi^(a). In many cases it is sufficient to only consider Pi^(a) since the 
other propagators Piy{<y) and Piz{<y) are related by a similarity transformation. In general, 
the method will depend on the error model under consideration; here we explicitly show 
how to construct this term in the amplitude and pulse-length error models and also in the 
addressing error model. 

Addressing errors: We wish to perform the evolution Pij.{a) using a product of imperfect 
square-pulse propagators M{9, 0). Recall that in this model rotations on the addressed qubit 
are error free MiO^cf)) = R{9,(j)), whereas on the unaddressed qubit the applied unitary 
depends on the systematic addressing error M{9,(f)) = R{9eN,4>)- Similarly, the sequence 
implementing Pix{(y) must resolve to the identity on the addressed qubit, while on the 
unaddressed qubit exp(— iaejv-f^x) + C){^'n) is applied. This behavior may be achieved by 
using a pulse sequence to implement Pi^- Let 

Sx{a) = M(27ra, -0„)M(27ra, 0«), (36) 

where a = [|a|/47r] is a integer number of 27r rotations and 0q = arccos(a/47ra). On the 
addressed qubit, Sx{a) resolves to the identity, whereas for the unaddressed spins, 

Sx{a) = M{2na, —(f)a)M{2na, (pa) = R{2TTaeM, —(f>a)R{,'^Ti^ci^N, 4>a) = Pix{c(), 
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Figure 6: Generation of the pure error term in the Solovay-Kitaev method, a) From Eq. (34), 
-f2z(tt/3) may be produced by the sequence Sx{—a)Sy{—(3)Sx{Q:)Sy{f3). On the Lie algebra 
the sequence corresponds to a closed rectangular path with an enclosed area of e^a(3. b) 
Alternatively, the rhombus construction may be used to generate P2z{ci) using four pulses. 



is applied. In the Lie algebraic picture, Sx{a) is composed of two vectors constructed so 
that their vector sum is —iaeNHx- Similarly, Sy{/3) = M{27ib, 7i/2 — (f)p)M{2rch, 7r/2 -|- (pj^) = 
Piy{f3). From these basic sequences, we may construct P2z{o:f3) using by using the balanced 
group commutator Eq. (34). In figure 6a we plot the sequence P2z{c((3) as a vector path on 
the Lie algebra. The sequence encloses a signed area e%a(3, which is denoted by the shaded 
rectangular figure. By tuning the rotation angles a and (3, one may generate a term which 
encloses any desired area, thus allowing the synthesis of an arbitrary pure-error term. 

At this point, we discuss a subtle feature of the addressing error model which at first 
appears to complicate the application of the SK method. Observe that on the unaddressed 
spin, the imperfect propagators may only apply small rotations (i.e. rotations by angles 
9eN)- If we restrict ourselves to sequences composed of resonant square pulses, then the term 
proportional to Piz{<y) may not be produced; we may not prepare such a term by similarity 
transformation (e.g. R{7i/2, 0)Sy{a)K^{7!'/2, 0)) since such an operation would either require 
a large rotation or if instead the transformation was carried out on the individual sequence 
propagators, the rotation axes would be lifted out of the H^-Hy plane. Similar arguments 
show that the Euler decomposition is also unavailable. 

Fortunately, this restriction is not as serious as it first appears; the SK method may be 
used provided that the sequence terms are chosen with care. We are ultimately saved by the 
orientation of the error terms in the Lie algebra. Using the BCH formula it is straightforward 
to show that for sequences composed of resonant pulses, the even-order error terms are always 
aligned along the —iH^ axis, whereas the odd-order terms are confined to Hy^ — Hy plane. 
Likewise, using only Sx{(y) and Sy{P), it is possible to generate correction terms that follow 
the same pattern. As a consequence, in this case it is possible to generate the correction 
terms f/rexp(— e^"'"^r2„_|_i)f/|, by carefully choosing the rotation angles and phases in the 
correction sequence -B„+i. 

As a instructive example, we shall derive a second-order passband sequence using the 
Solovay-Kitaev method. We begin by calculating the Magnus expansion for the seed se- 
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quence Mski(6',0) = [/Texp(e^r22 + e%fl3 + ■ ■ ■ ) where the target operation Ut = t for 
the unaddressed qubit. To cancel the second order term, we simply need to apply the 
inverse of exp(e^r22) = exp(— i27r^e^ sin(20sKi)-f^z)- The planar rotation T = 1, since 
B2 = P2z(— 27r^ sin(20sKi)) is already oriented in the correct direction. One possible choice 
for B2 is the sequence 

B2 = Sa;{-27icos(f)sKi)Sy{27ism(j)sKi)Sx{2iTcos(f)sKi)Sy{-27!-sm(j)sKi)- (37) 

The sequence Msk2{0, 0) = B2Mski{0, 0) = Ut + 0{e%) corrects addressing errors to second 
order. We denote an nth order compensating sequence produced by the Solovay-Kitaev 
method using SKI as an initial seed as SKn; here we have produced an SK2 sequence. 

More efficient constructions for the correction sequence B2 exist. Observe that we may 
directly create the pure error term P2z{l) by using four pulses in the balanced group 
commutator arrangement P2z{l) = M"'"(27rc, 0^)M"'"(27rc, 0l^)M(27rc, 0^)M(27rc, 01^), where 
c = [|7|/47r] and the phases are chosen to be 0^ = arcsin(7/47r^c^) and cp'^ = (1 — sign7)7r/2. 
The phase (p' is only necessary to ensure that the construction also works for negative 7. 
We call this arrangement the rhombus construction. In figure 6b we plot this sequence as a 
vector path on su(2). In this construction the magnitude of the error term is tuned by ad- 
justing the phase (p^, i.e. adjusting the area enclosed by the rhomboidal path of the sequence 
in the Lie algebra. The rhombus construction has the advantage of requiring half as many 
pulses as the standard method. We will now use this construction to produce an alternative 
form for SK2. Let 

B'^ = M(27r, 0^ + 7i)M{2n, 0)M{2n, 0^)M(27r, vr) (38) 

where (p^ = arcsin(sin(20sKi)/2). Then M^^2{e,0) = 5^Mski(^,0) = Ut + 0{e%). 

Amplitude / pulse-length errors: We now turn our attention to the compensation of 
amplitude and pulse-length errors. When considering these error models, the imperfect 
propagators of the form M{9,(p) = R{9{1 + eA),<P)- In this case we may also construct 
Pix using Eq. (36) since the imperfect 27r rotations reduce to M(27ra, 0) = R{2TTaeA, (p) 
and Sx{oi) = Pix{(y). As a consequence, if the initial seed sequence Wn is a passband 
sequence then Wn+i is also a passband sequence. We note however, that in this error model 
one has more fiexibility in the synthesis of the correction sequences B^+i; since now the 
imperfect propagators apply large rotations, then we may perform similarity transformations 
of sequences by directly implementing the required pulses. Notably, we may use imperfect 
propagators of order 0{e^) to apply a desired transformation at a cost of an error of order 
0{el+^), for example M{e,(p)Pj,M\e,<p) = R{e,<p)Pj,R\e,(p) + 0{e\). 

The SK method may be used to calculate higher order compensation sequences for the 
amplitude error model. We begin by calculating the relevant Magnus expansion for the seed 
sequence Mski(^, 0) = M(27r, -0sKi)M(27r, 0ski)M(^, 0) = UTexp{e\n2 + e\n3 + - ■ ■ ) where 
now in the amplitude error model the target operation Ut = R{6,0). To cancel the second 
order term we must apply the inverse of Ut ewi^A^^W^, = exp(— i27r^e^ sin(20sKi)-ffz). This 
is precisely the same term that arose previously for addressing errors and the error may be 
compensated the same way. In fact every nth order SK sequence is passband and works for 
both error models. 
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4.2 Wimperis / Trotter-Suzuki sequences 

In this section we study a second family of fully compensating pulse sequences, first discov- 
ered and applied by Wimperis [52], which may be used to to correct pulse length, amplitude, 
and addressing errors to second order. These sequences have been remarkably successful and 
have found extensive use in NMR and quantum information [54, 55, 56]. Furthermore, the 
Wimperis sequences may be generalized to a family of arbitrarily high order sequences by 
connecting them to Trotter-Suzuki formulas [51]. An analogous composite sequence com- 
posed of rotations in SU(4) may be used to correct two-qubit operations [57, 58]. We study 
the two-qubit case in section 5.1. 

Narrowband behavior: We begin with the problem of identifying narrowband sequences 
which correct addressing errors. In the addressing error model operations on addressed 
qubits are error free, whereas on the unaddressed qubits the imperfect propagators takes the 
form V{u(t)) = U{eNu(t)) and Ut = 1- Further, we shall constrain ourselves to sequences 
composed of resonant square-pulse propagators. Specifically, we search for arrangements of 
four pulses that eliminate both the first and second-order Magnus expansion term for the 
imperfect propagator. 

Before explicitly describing the construction of the Wimperis sequences, we digress shortly 
to point out a certain symmetry property which may be used to ensure that e1fQ2 = 0. 
Consider the group product of propagators of the form, 

UieNu{t)) = U2U3U2U1, Uk = exp{-ieNtkUk ■ H), (39) 

with the added condition that the first-order expansion term for U{eNu{t)) has been already 
eliminated, i.e. eN^i = —i^Nitsu^ + 2t2'"2 + tiUi) ■ i? = 0. In this arrangement, the second 
and fourth propagators are identical; this pulse symmetry along with ti^i + 2t2ii2 + ^3'"3 = 
eliminates the second-order term. 



e 



2 4 j 
^^2 = ^^^[-it.n,-/f,-it,ii,-if]=0. (40) 



As a consequence, U{eNu{t)) = 1 -|- 0{e%). Alternatively, one may regard the product 
T = U2U3U2 as a second-order symmetric Trotter-Suzuki formula for the inverse operation 
f/| [47] which approximately cancels the undesired rotation Ui. Considering the control fields 
applied during the application of the corrector sequence T, the applied control Hamiltonian is 
symmetric with respect to time inversion. By a well known theorem, all even-order expansion 
terms produced by a time-symmetric Hamiltonian cancel [43, 59] (i.e. ej^^2j = 0, for all 
positive integers j). Thus by implementing symmetric corrector sequences T2J, it is sufficient 
to only consider the cancellation of the remaining odd-order error terms. In section 4.2.1 
we will inductively develop a series of corrector sequences of increasing accuracy based on 
symmetric Trotter-Suzuki formulas [47, 51]. 

The cancellation of the second-order term may also be inferred from geometric consid- 
erations on the Lie algebra. To be concrete, consider a sequence of the form Eq. (39) 



4 COMPOSITE PULSE SEQ UENCES OiV S U (2) 



27 





Figure 7: a) Vector path followed by N2 on the Lie algebra, b) Trajectory of an unaddressed 
spin during an N2 sequence, using imperfect rotations of the form M{6,(f)) = R{6eA,4')- c) 
B2 correcting an amplitude error, using imperfect rotations of the form M{6, 0) = R{6{1 + 
e^), 0)- In these plots e^ = ^a = 0.2. 



where 



e^tiUi ■ H 
eNtsUs ■ H 



7reiv(cos 0N2-f^x + sin (pmHy) 
2tc eN (cos (f)mH^ - sin0N2-f^y), 



(41) 



and the phase is 0n2 = arccos(— ^/47r) is chosen so that tiUi + 2^21*2 + taU^ = 0, i.e., the 
vectors form a closed path on the dynamical Lie algebra. Figure 7a is a diagram of these 
vectors on su(2). In su(2), the Lie bracket is equivalent to the vector cross product [35], 
therefore we may interpret ||e^n2||HS as the signed area enclosed by the vector path on the 
Lie algebra. We note that the sequence under study encloses two regions of equal area and 
opposite sign, ensuring that the second order term is eliminated. 

With this insight, the construction of a second-order narrowband sequence is straightfor- 
ward. The vectors Eq. (41) correspond to the sequence 



MmiO, 0) = M(7r, 0N2)M(27r, -0N2)M(7r, 0n2)M(^, 0). 



(42) 



Wimperis refers to this sequence as NBl, and indeed, this is the established name in 
the literature. In the current article, we label this sequence N2 in anticipation of the 
generalization of this form to N2j, which compensates addressing errors to 0{2j). We 
use this language to avoid confusion with other established sequences, namely NB2, NB3, 
etc [52]. The N2 sequence may be used to compensate addressing errors. For unad- 
dressed qubits M{9,(j)) = R{9eN,4>) a^id thus from Eq. (40) and Eq. (41) it follows that 
Mn2(6',0) = 1 + 0{e%). Thus on unaddressed qubits the sequence performs the identity 
operation up to second-order. Furthermore, for addressed qubits M{9, 0) = R{9, 0) and 
therefore Mn2(6', 0) = R{9, 0). As a result, when the sequence Mn2(^, 0) is used in the place 
of the imperfect operation M{9,0), the discrimination between addressed and unaddressed 
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spins is enhanced. In figure 7b we plot the magnetization trajectory for an unaddressed 
qubit under an N2 sequence. 

Broadband behavior: As previously discussed, broadband sequences are best suited for 
correcting amplitude or pulse-length errors. In the following we will explicitly consider the 
amplitude error model, where M{9,(j)) = R{6{1 + eyi),0) = R{9,(f))R{9eA,(p)- Although a 
pulse sequence may be studied by considering the interaction frame propagator as described 
in section 3.2.1, the method originally used by Wimperis is simpler. Wimperis' insight was 
that for lowest orders the toggled frame could be derived geometrically, using the relation 
R{9, —(f)) = R{7i,(j))R{9,3(f))R{TT,(f) + tt). Consider the apphcation of the target gate Ut = 
R{9, 0) using the sequence 

Mb2{9, 0) = M(7r, 0B2)M(27r, 30B2)M(7r, 0b2)M(^, 0) (43) 

where 032 = arccos(— 6'/47r). This is the Wimperis broadband sequence, traditionally called 
BBl but here denoted as B2. When rewritten in terms of proper rotations one obtains, 

Mb2(^, 0) = R{7ieA, 0B2) [R{71, 0B2)i?(27reA, <pB2)RH7r, 0B2))i?(vreA, MRid^A, 0)R{9, 0) 
= R{7ieA, 0B2)^(27re^, -0B2)^(vreA, MRi^^A, 0)] R{9, 0) 

where the identity R{27i, 30b2) = —1 = R{t^, 0B2 + vr)i?(7r, 0b2 + ti") was used. Let Q denote 
the quantity enclosed in square brackets, so that we may write Mb2{9,0) = QUt- Observe 
that Q is precisely the form considered previously in N2. From our previous result we may 
conclude Q = 1 + 0{e\) and Mb2(^,0) = R{9,^) + 0{e\). As a result, when Mb2(^,0) is 
used in the place of the imperfect operation M{9, 0) the effect of the systematic amplitude 
error is reduced. Note that errors for the B2 and N2 sequences follow equivalent paths on 
the Lie algebra in their respective interaction frames. In figure 7c we plot the magnetization 
trajectory for a qubit under a B2 sequence with amplitude errors. 

Passband behavior: In some cases, it is convenient to have a passband pulse sequence 
that corrects for both addressing errors and for amplitude errors, as we saw with the 
Solovay-Kitaev sequences. The passband Wimperis sequence P2 is simply two Solovay- 
Kitaev correction sequences in a row where the order of pulses is switched, Mp2(^,0) = 
M(27r,0p2)M(27r, -0p2)M(27r,-0p2)M(27r,0p2)M(^,O), with 0P2 = arccos(-^/87r). The 
switching of pulse order naturally removes the second order error term. One may verify that 
this sequence works for both addressing and amplitude errors. 

4.2.1 Arbitrarily accurate Trotter-Suzuki sequences 

The Wimperis sequences rely on a certain symmetrical ordering of pulse sequence propagators 
to ensure that even order Magnus expansion terms are eliminated. We may further improve 
the performance of these sequences by taking advantage of additional symmetries that cancel 
higher-order terms. In the following, we shall show how by using symmetric Trotter-Suzuki 
formulas [46, 60, 47] a family of arbitrarily accurate composite pulse sequences may be 
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constructed [51, 53]. The Lie algebraic picture along with the Magnus and BCH series will 
be important tools in this process. 

Symmetrized Suzuki formulas: Before discussing the particular form of these sequences, 
it is helpful to briefly mention a few important results regarding symmetric products of time- 
independent propagators [47] . Given a series of skew-Hermitian time- independent Hamilto- 
nians {Hi, H2, . . . , Hm} such that ^™^ Hi = Ht, the BCH expansion tells us 

m / 00 \ 

J] exp{XHi) = exp A^t + J^ A"fi„ , (44) 

1=1 \ n=2 / 

where A is a real parameter and the expansion terms Qn depend on the specific ordering of 
the sequence. If we choose to apply the operators in a time-symmetric manner such that 
exp(Aifj) = exp{\Hm+i_i), then the symmetry of the pulse removes all even-order terms. 
For symmetric products, we have 

Yl exp(XHi) = exp [ XHt + ^ A^^-^^sj-i ) • (45) 

symmetric \ j=2 / 

An important observation concerning the elimination of the remaining odd terms was 
made by Suzuki [47]. Provided that Hi = PiHr + {'PiY^~^HB^ where the coefficients Pi 
are real numbers, there exist certain choices of coefficients such that YIT ^i- ~ ^t- This 
requires that YIT Pi ~ ^ ^^^ STfe)^"'"^ ~ 0- "^^^ situation is considerably simplified if 
we restrict ourselves to sequences composed of just two kinds of propagators, f/i = exp(ifi) 
and U2 = exp{H2). In this case the previous expression simplifies to nipi + ^2^2 = 1 and 
iT'iPi' + ^2P2~ = 0, where the integers ni and n2 are the number of Ui and U2 pulses 
required to produce a sequence that is independent of Hb up to 0{p^^^^). We solve for a set of 
coefficients by setting p2 = — 2j9i and n2 = 1, thus yielding ni = 2^-'"^ and pi = 1/(2^-'"^ — 2). 

Combining these observations, if W2j-2{p) is a (2j — 2)th approximation of exp{pHT) and 
Q2J-1 = P^-'^^Hb where Hb is independent of p, then we can construct W2j{l) from the lower 
approximations, W2j{l) = (W^2j-2(p))^ ^ W^2j-2(— 2p)(VF2i-2(p))^ ^ • As a result, Suzuki 
formulas provide a path of producing higher-order sequences from a symmetric combination 
of lower-order sequences. Notice that the symmetric decomposition is used to keep the 
even-order terms zero. 

Passhand behavior: In the following we will seek to generalize the second-order passband 
sequence P2 to an arbitrarily accurate passband sequence P2j. Our goal is to develop a 
correction sequence T2j{k,(f)) = exp(i6'e^ifx) + ^{^a^ ) that cancels the unwanted rotation 
of a M{9, 0) operation. 

When considering passband sequences, it is convenient to use rotation angles that are 
integer multiples of 2tt. Let us define the triangular motif, 

Ti(A;,0) = M(2A;7r, -0)M(2A;7r, 0) 

= exp(— i4A;7reyiCos0ifx — i(27rA;eA)^cos</)sin0ifz) + O(eJ^). (46) 
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Figure 8: Vector path followed by P4 on the Lie algebra. 



Observe the passband sequence SKI may now be written as Mski{0, 0) = Ti(l, (f)sKi)M{9, 0) 
and where the phase 0ski is as defined previously. The remaining second-order term is odd 
with respect to (p, and is related to the vector cross product on the Lie algebra. A second- 
order sequence may be constructed by combining two Ti{k,(j)) terms so that the correction 
sequence is symmetric and the cross product cancels. Let us define the symmetrized product 



T2(fc,0) = Ti(A;,-0)Ti(fc,0) 



(47) 



where the length p = — {Skn/ 9) cos cf), the target Hamiltonian Ht = iOeAHx, and Hb = 
e^Q^/p'^ is the remaining term we wish to cancel. ^^3 is a function that depends on k and 
such that Hb is a function of but not of the length scale k. For fixed and variable k, this 
makes Hb independent of p. Observe that the passband sequence P2 may now be written 
as Mp2{9,0) = T2(l, 0P2)M(6', 0), where again 0P2 = arccos(— ^/Svr). 

Our strategy is to use a Suzuki formula to construct higher-order T^ik^cj)) (that is, 
T2j{k,(j)) for j = 2) using a symmetric combination of (ni = 8)-many exp{Hi) = T2{k,(f)) 
sequences and a single exp(if2) = T2{—2k,(f)) sequence, yielding 

(48) 



T4(A;,0)= (T2ik,(l))) T2(-2fc,0)(^T2(fc,0) 

In order to produce to required correction term, the parameters {k, 0) must be chosen such 
that {rii — 2)p = 1. The fourth-order passband sequence P4 is Mp4^{9, 0) = T4(l, 0P4)M(6', 0), 
where 0P4 = arccos(— 6*74871). In figure 8 we plot the vector path followed by P4 on the Lie 
algebra. This result may be further generalized. To produce T2j(/c,0) requires {rii = 2^-^"^)- 
many T2j-2{k,(j)) sequences and a single T2j_2(— 2/c, 0) sequence in the symmetric ordering, 

(49) 



T2,(A;,0) = (t2,_2(/c,0)) T2,_2(-2A;,0)(t2,„2(A:,0) 

We then fix so that the first-order term cancels the unwanted rotation, yielding 

9 



>'P2j 



arccos 



27r/, 



(50) 
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Figure 9: Infidelity of the Trotter-Suzuki sequences B2j, P2j, and N2j. In the amplitude 
error model (top row) M{9,(f)) = M{9{1 + eA),4>) and Ut = -R(vr/2,0). In the addressing 
error model (bottom row) on the unaddressed qubits M{9,(j)) = M{6eN,<t') aiid Ut = 1, 
while on the addressed spins i?(7r/2, 0) is applied. Each error model establishes a separate 
preferred interaction frame; when transformed into the appropriate pictures, the B2j and 
N2j sequences are homologous. The passband sequences P2j can correct both amplitude 
and addressing errors at the cost of reduced efficacy. 



where fj = (2^-' ^ — l)/j-i and for the sequence P2j, /i = 4. Then the 2jth-order passband 
sequence P2j is Mp2j(^,0) = T2j{l,(j)P2j)M{e,0). 

The same method can be used to develop exclusively narrowband or broadband sequences 
called N2j and B2j respectively [51]. This requires redefining the bottom recursion layer 
T2{k, (f)) to have either narrowband or broadband properties. For N2j, f\ = 2 and T2{k, 0) = 
Ti{k/2,(j))Ti{k/2,—(j)). For B2j /i = 2, but T2 is slightly more complicated; when k is 
even T2(fc,0) = Ti(/c/2, — 0)Ti(fc/2, 0) just like N2j, however when k is odd and T2(fc,0) = 
M {kn , (p) M {kir , 3(j)) M {kir , 3(j)) . In figure 9, we compare several of the generalized Trotter- 
Suzuki sequences to the ideal unitaries Ut = R{n/2, 0) in the case of amplitude errors (top 
row) and Ut = t in the case of addressing errors on unaddressed qubits (bottom row). 
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4.3 CORPSE 

So far, the sequences considered here have been designed to correct systematic amphtude and 
addressing errors. The correction of errors arising from an inaccurate tuning of the control 
field are also of practical interest. The treatment of detuning errors is similar in principle to 
the error models already considered, however in practice the construction of compensating 
sequences is complicated by the noncommutivity of the ideal Hamiltonian and the erroneous 
Hamiltonian generated by the control distortion [u{t) ■ H, Su{t) ■ H] ^ 0. 

Fully compensating pulse sequences for detuning errors were originally studied by Tycko 
[61] and later generalized by Cummins and Jones into the humorously named ROTTEN 
(resonance offset tailoring to enhance nuition) [62] and CORPSE (compensating for off- 
resonance with a pulse sequence) [63, 64] family of sequences. Cummins and Jones initially 
derived the sequence using quaternion algebra to represent simple rotations, and optimized 
the angles to eliminate the first order effects of the detuning error. The CORPSE family 
of sequences has found application in NMR [63] and SQUID [65] experiments. Here we 
reexamine CORPSE using the techniques outlined in section 3.2.1. 

CORPSE is a sequence which performs a compensated rotation about the H^ axis. Fol- 
lowing Cummins and Jones, the sequence is constructed from three square pulses, which in 
the case of perfect resonance, induce rotations about the H^, —Hy^, and H^ axes, sequentially. 
The pulse sequence is parametrized by the piecewise constant control function 

{Mx < t < ti 

- Mx ti<t<t2 , (51) 

Mx t2 <t <T 

where u^ is a constant H^ control amplitude and the tj are times at which the field direction 
is switched. This construction is particularly amenable to analytic methods, as the ideal 
control Hamiltonian H{t) = Ux{t)H^ commutes with itself at all times. In the absence 
of systematic detuning errors the control Hamiltonian produces the following ideal unitary 
evolution, 

f/(u,(t')x; t, 0) = exp(-i^(t)i7x) = i?(^(t), 0), ^{t) = f dt'u^{t'). (52) 

Jo 

Over the entire time interval < t < r, the pulse sequence produces the gate Ut = R{0, 0), 
where 6 = '&{t). 

In the presence of an unknown detuning error du = e^z and the rotation axis is lifted in 
the direction of the H^ axis on the Bloch sphere. Recalling section 3.2.1, a Magnus expansion 
may be used for the interaction frame propagator ^/^(e/jz; r, 0), which produces the evolution 
generated by the systematic detuning error. Combining Eq. (17) and Eq. (52) the first-order 
term is 

eD^i{r,0) = -ieD I dtcos{-d{t))H^ + sm{'d{t))Hy. (53) 

Jo 
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Direct integration yields, 

e^fii(r, 0) = -'-^[sm{e,)H, + (1 - cos{ei))Hy] + ■■■ 

+ ^[(sin(^i - 62) - sm{ei))H, + (cos(^i) - cos(^i - 02))^y] + ■ ■ ■ (54) 

-^[(sin(^i -e, + ^3)- sin(^i - 02))H,+ (005(^1 - ^2)- cos(0i - ^2 + ^3))i^y], 

where 6^ = ^{tk) — i^ltk-i) are the effective rotation angles applies during the kth square 
pulse. At this point, we may interpret each of the terms of Eq. (54) as vectors on the 
dynamical Lie algebra. Figure 10a is a diagram of these vectors on su(2). In order to 
eliminate the first-order expansion term, the rotation angles 6k must be chosen so that the 
vectors must form a closed path. 

Then by choosing the rotation angles to be 



9, -- 


= 27ini + 9/2 - arcsin (sin(^/2)/2) 


02 -- 


= 271722 - 2 arcsin (sin(^/2)/2) 


03 - 


= 27^3 + 6/2 - arcsin (sin(e/2)/2) 



(55) 

where rii, 77-2 and n^ are integers, we find that both eDQi{r,0) = and 9 = {9i — 62 + 6*3 
mod 27r). The extra factors of 27r are added so that the individual pulse rotation angles may 
be made positive. In principle, any choice of integers is sufficient to compose a first-order 
sequence, however the choice ni = n2 = 1 and 113 = minimizes the remaining second-order 
term while still producing a positive set of rotation angles [64]. The CORPSE family of 
sequences, 

V{uit)) = exp(-i(^3i^x + ezp03^z)) exp(-i(-02i^x + eD^2i^z)) exp{-i{9iH, + tD9iH,)) 
Mci(0,O) = Af(^3,O)Af(e2,7r)M(0i,O)=i?(^,O) + O(e^) (56) 

are fully compensating first-order sequences. In the presence of an unknown detuning error, 
a CORPSE sequence, Mci(6',0), may be implemented in the place of the simple rotation 
R{9, 0) and the erroneous evolution is suppressed. In figure 10b, we plot the magnetization 
trajectory for a qubit under a CORPSE sequence with detuning error eo = 0.2. 

4.3.1 Arbitrarily accurate CORPSE 

We turn our attention to sequences which compensate detuning errors to arbitrarily high 
order. Once again, the Solovay-Kitaev method may be used to construct higher-order se- 
quences, now using CORPSE as the seed sequence. This problem was first studied by Alway 
and Jones [53]. Recall that in the Solovay-Kitaev method, one synthesizes a correction se- 
quence Bn+i in two steps. First, one generates a propagator P{n+i)z{i) with an amplitude 
proportional to the leading-order error term. Second, the similarity transformation f/j'T"'' is 
used so that -B„+i = UtT'^ P{n+i)z{i)'^U}p cancels the leading error term of the seed sequence 
(see Eq. (35)). Eq. (34) shows that the term P(n+i)z{,C) may be constructed recursively 
using a product of first-order Pi's. Therefore the Solvay-Kitaev method may be extended to 
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Figure 10: a) Vector path followed by CORPSE on the Lie algebra, with the choice of 
parameters rii = n2 = n^ = 0. Each vector g^^ corresponds to a term in Eq. (54). b) 
Trajectory of a spin under a CORPSE sequence for Ut = -R(7r/2, 0), with rii = n2 = I and 
ris = chosen to produce positive angles. The sequence is constructed of imperfect rotations 
of the form M(^fe,0) = exp(-i(0fc/7x + toOkH^)), with en = 0.2. 



correct detuning errors if: (1) procedures for generating Pi^, Piy, and Piz using imperfect 
pulses have been found, and (2) when a method for applying the similarity transformation 
UtT^ using imperfect pulses has been identified. This task is complicated by the difficulty 
of generating inverse operations for imperfect pulses affected by detuning errors. 

Earlier, the first-order Pi's were applied using a composite pulse sequence (see Eq. (37)). 
We employ a similar strategy here. Let 



Sz{a) = M(a/2, + 7r)M(a/2, 0) = exp{-iaeH^) = Pi^(a). 



(57) 



Then we may apply ^12(0;) by implementing the sequence Sz{a) instead. What remains is to 
develop sequences that apply Pix{a) and Piy{a). Observe that if it were possible to perform 
ideal rotations (where R\0,(l)) = R{6,(j) + vr)), then one could easily implement Sx{oi) = 
R{TT/2,TT/2)Sz{a)R^ (71/2,71/2) and Sy{a) = R^{7r/2,0)Sz{a)R{7i/2,0) in the place of Pi^.(a) 
and Piy{a). However in the presence of systematic detuning errors, we may not apply 
this transformation using a simple imperfect pulse since M'^{6,(j)) 7^ M{6,(f) + vr). We may 
avoid this complication by using a CORPSE sequence to approximate the transformation to 
sufficient accuracy. We use the first-order CORPSE sequence, Mci(6',0) to implement the 
target operation Ut = R{G, 0)- Then we may write 



Sxia) = Mci(7r/2, 7r/2)S,(a)Mci(7r/2, 37r/2) = exp{-iaeH^) + 0{e^) = P^a) (5^ 



and 



Sy{a) 



Mci(7r/2,7r)5,(a)Mci(7r/2,0) = exp(-iae/7y) + 0{e^) = P, 



ly 



a 



(59) 



Furthermore, by careful choice of 6 and 0, we can implement Pi^(a) for any axis rj about any 
angle a. Following the Solovay-Kitaev method (Sec. 4.1.1), we can then construct PjziO to 
any order j. 
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What remains is to implement the similarity transformation by f/j-T'''. However we note 
that the available approximate rotations Mci(^,(/>) = R{9,(j)) + 0{e£)) are only accurate 
to first order. Hence, if we were to attempt to apply a similarity transformation on a 
higher-order PjziO term using CORPSE sequences, additional second-order errors would 
be introduced. This difficulty may be avoided by applying the appropriate rotation at the 
level of the first-order Pi operations. As an instructive example, we work through the 
process for the second-order correction. The correction sequences may be constructed in 
a rotated coordinate system, determined by the basis transformation H' = Ut'^'^ H ^TUj. 
for /i G {x, y, z}. One may calculate four rotations R{9±^i,(f)±^i)^ that map the H^ axis 
to the if±x' and H±yi axes H±^j_i = R{0±^j,i,(t)±fj,')H.^R\6±^i,(l)±^i) and also Pi±^i{q.) = 
R{d±fj.',4>±iJ.')Piz{c()R'^{0±ij.',4>±fi')- This transformation may also be applied using a first- 
order CORPSE sequence, since the resulting error is absorbed into the second-order error 
term in Pi±^'{<y)- Once transformed terms have been obtained, the correction sequence may 
be constructed in the usual manner 

B2 = P2Ai) = PiA-Vl)Piy'{-Vl)PiAVl)Piy'{Vl) 

xMci{e,>, 0x0 A^( V^)^ci(^x', 0x' + AMci{9y,, (i)y')Pi,{^i)Mci{ey,, (Py, + 7r). 

We can then define a second order CORPSE sequence as Mc2(^, 0)= BiMciiO, 0)= R{e, 0) + 
0(e|)). One can then continue using Solvay-Kitaev techniques to remove the detuning error 
to all orders. 

4.3.2 Concatenated CORPSE: correcting simultaneous errors 

So far, we have considered the problem of quantum control in the presence of a single 
systematic error, while in a real experiment several independent systematic errors may affect 
the controls. In many situations one error dominates the imperfect evolution; it is appropriate 
in these cases to use a compensating sequence to suppress the dominant error. However it is 
also important to consider whether a sequence reduces the sensitivity to one type of error at 
the cost of increased sensitivity to other types of errors [64] . Such a situation may occur if 
an error model couples two error sources. Consider a set of control functions {^/^(t)} where 
each control is deformed by two independent systematic errors. A general two-parameter 
error model for the control v^{t) = /^[■u(t); e^, ej] may be formally expanded as 

f^[u{t)-e„e,] = /^[t,(t);0,0] + e,— ^[ii(t);0,0]+e,— ^[ii(t);0,0] + 

rj2 

2e,ej-^^^Uuit); 0, 0] + Oie^ + e]). (60) 

In practice, we need only concern ourselves with the first few expansion terms, since in most 
physically relevant error models the higher order derivatives are identically zero. Following 
the same reasoning employed in section 2.1, we may decompose the imperfect propagator as 
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V{uit)) = U{u{t))U\v{t) — uit)), where U{u{t)) is the ideal propagation in the absence 
of errors and the interaction frame propagator U\v{t) — u{t)) represents the evolution in- 
duced by the systematic errors. Formally U\v{t) — u{t)) may be studied using a Magnus 
expansion, however this method is usually impeded by the complexity of the Magnus series. 
The determination of sequences which compensate simultaneous errors is currently an unre- 
solved problem, although some progress has been made by considering concatenated pulse 
sequences [31]. As an example relevant to an NMR quantum computer, we now study error 
models where two systematic errors occur simultaneously and show that these errors can be 
compensated by concatenation of pulse sequences. 

Simultaneous amplitude and detuning errors: When studying the control of qubits based 
on coherent spectroscopy methods, it is natural to consider situations where systematic 
errors in the field amplitude and tuning are simultaneously present. From the NMR control 
Hamiltonian Eq. (9) it is straightforward to derive the error model, fx/y = Uy^/y{t){l + eA) and 
Vz = Uz(t) + en- In this case, the independent amplitude and detuning errors are decoupled 
and the parameters e^ and eo do not affect the same control. 

Simultaneous pulse-length and detuning errors: Likewise, simultaneous errors in the pulse 
length and field tuning may occur. Again, it is simple to derive the joint error model, 
v^/y = Mx/y(i)(l + ct) and v^ = {u^if) + e£,){l + ej-). Unlike the previous case, this error 
model couples the parameters ey and e^, since they both effect the H^ control function. 

It is sometimes possible to produce a sequence that compensates two simultaneous errors 
by pulse sequence concatenation. Recall from section 4.2 that the Wimperis sequence B2 is 
a second-order compensation sequence for both amplitude and pulse-length errors. However 
if only detuning errors are present, then the response of a B2 sequence to pure detuning 
errors is similar to that of a native pulse Similarly, a CORPSE sequence would eliminate 
the first order term produced by the detuning offsets but not correct amplitude errors. 
The sequence B2C0RPSE corrects both errors independently by concatenating the B2 and 
CORPSE sequences. B2C0RPSE is composed of three B2 subsequences that form a larger 
CORPSE sequence 

Mb2Ci(^,0) =MB2(^3,0)MB2(^2,Vr)MB2(^l,0), (61) 

where the angles {^1,^21^3} are given in Eq. (55). At the lower level of concatenation, 
B2 sequences are used to synthesize rotations robust to pure amplitude or pulse-length er- 
rors, whereas on the higher level the CORPSE construction compensates pure detuning 
errors. Figures lla-d are plots of the infidelity of plain pulses, and the CORPSE, B2, and 
B2C0RPSE sequences in the presence of simultaneous amplitude and detuning errors. Sim- 
ilarly, figures lle-h are infidelity plots for simultaneous pulse- length and detuning errors. As 
expected, the CORPSE sequences improve the accuracy of gates with respect to detuning 
errors, but offers little improvement against either amplitude or pulse length errors. Unsur- 
prisingly, this behavior is inverted for the B2 sequences. The B2C0RPSE sequence performs 
well for either amplitude/pulse-length or detuning errors; in the presence of simultaneous 
errors the performance diminishes, yet the fidelity of the applied gate is still vastly improved 
over uncompensated pulses. 
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Figure 11: Infidelity of a,e) Plain pulses, b,f) CORPSE c,g) B2, and d,li) a concatenated 
B2C0RPSE sequence in the presence of simultaneous amplitude and detuning (left column) 
and pulse-length and detuning (right column) errors. The target rotation is Ut = -R(7r,0). 
The dashed contour corresponds to an infidelity of 0.01, while the remaining contours are 
plotted at 10% intervals. 
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4.4 Shaped pulse sequences 

Thus far, we have considered sequences composed control functions M^(t) which are piecewise 
constant over the pulse interval. This construction is particularly convenient from a sequence 
design perspective, although in practice instrumental shortcomings will frequently distort the 
pulse profile. Also, in some applications, the rectangular profile is non-ideal. For example, in 
the frequency domain the square pulse corresponds to a sine function, whose local maxima 
may complicate the control of certain systems. Finally, we note that in some systems, such 
as in superconducting Josephson junction qubits, bandwidth requirements forbid the sudden 
switching of control fields (i.e., place a limit on the derivative Im^I < v}^^^). In these cases, it 
is desirable to consider shaped pulse sequences composed by a set of continuous differentiable 
control functions [66]. In the present article, we shall only consider shaped pulse sequences 
which are also fully-compensating (class A), and therefore appropriate for use in a quantum 
processor. Specifically, we study shaped pulse sequences which also compensate systematic 
errors. We shall see that the additional fiexibility admitted by shaped pulses frequently 
produces superior sequences. 

Once given a pulse waveform, we may verify that it is compensating for a particular 
error model by computing the Magnus expansion in the appropriate interaction frame (see 
section 3.2.1); however, we emphasize these methods require the evaluation of successive 
nested integrals (e.g., Eq. (15)), and are difficult to analytically implement beyond the 
first few orders [67]. Furthermore, the inverse problem (solving for the control functions) is 
especially difficult except for in the most simple cases. For this reason, various numerical 
optimization methods have become popular, including gradient-accent techniques, optimal 
control methods [68, 69, 70], and simulated annealing [71, 72]. Methods which use elements 
of optimal-control theory merit special attention; in recent years the GRAPE [73, 74] and 
Krotov algorithms have been especially successful in pulse design, and has been appUed to 
NMR [68], trapped ions [75], and ESR [76]. The main advantage of the GRAPE algorithm 
is an efficient estimation of the gradient of the fidelity as function of the controls, which then 
enables optimization via a gradient-accent method. A recent review of these methods may 
be found in Refs. [77, 78]. 

To demonstrate the relative performance of shaped sequences, we consider the continuous 
analogs of the CORPSE sequence that modulate Ux{t) to compensate detuning errors. From 
Eq. (53) observe that the first-order compensation condition eDQi{T,0) = is met when 
J^ dtcos{^{t)) = J^ dtsm{{}{t)) = 0. What remains is to solve for a set of control functions 
Mx(t) which generate the target operation Ut while remaining robust to the detuning error. 
One popular scheme is to decompose the control functions as a Fourier series [72, 79] over 
the pulse interval, 

oo „ 

Ux{t) = uy^ancos incu it — —]] + bnSin (nu it — -]] , u = — . (62) 

n=0 

This has the advantage of specifying the controls using only a few expansion parameters {a„} 
and {bn}- Also, the series may be truncated to avoid high frequency control modulations that 
are incompatible with some control systems. Moreover, by choosing each of the expansion 
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Figure 12: a) Vector path on the Lie algebra traced out by the interaction frame Hamiltonian 
H^{t) = eciH^it) for the shaped sequence Si(7r/2). b) Magnetization trajectory for a qubit 
under an Si(7r/2) sequence for the target operation Ut = -R(7r/2,0), for the exceptionally 
large detuning error cd = 1. 



terms 5„ = 0, the resultant sequence may be made symmetric with respect to time reversal. 
We note that not all compensating sequences completely eliminate the first order error term; 
in some sequences the leading order errors are highly suppressed rather than completely 
eliminated. This behavior is more common in sequences obtained from numerical methods. 

Historically, among the first class A shaped pulses developed belonged to the U-BURP 
(band-selective, uniform response, pure-phase) family designed by Geen and Freeman [72, 
25], using simulated-annealing and gradient decent methods. Soon after, Abramovich and 
Vega attacked the same problem using approximation methods based on Floquet theory 
[80]. These sequences are designed to correct detuning errors over an extremely large range; 
however, they fail to completely eliminate the leading order error terms. In the design of 
compensating sequences for quantum computing, the emphasis has been on the synthesis 
of composite rotations of extraordinary accuracy over a narrow window of errors. More 
recently, Steffen and Koch considered shaped Gaussian modulated controls specially designed 
for superconducting qubit manipulations [81]. Pryadko and Sengupta also derived a family of 
sequences {Si{(f)o) ,S2{4>o) ,Qi{(po) ,Q2{(Pq)) using a semi-analytic method based on the Magnus 
expansion and average Hamiltonian theory [82, 67]. By design, these sequences eliminate 
the leading-order error terms in a manner similar to the square-pulse sequences previously 
discussed. 

It is interesting to compare the performance of shaped sequences to a sequence of square 
pulses, such as CORPSE. Many of these sequences compensate detuning errors by modu- 
lating the amplitude of the control function u^{t). In figure 13a pulse shapes are given for 
CORPSE, U-BURP, Si(7r/2), S2(7r/2), Qi(7r/2), and Q2(vr/2) sequences that implement the 
target operation Ut = -R(vr/2,0). The relative performance of these sequence over a wide 
range of systematic detuning error is given in figure 13b. Clearly, shaped pulses outperform 
CORPSE over a wide range of field tunings; however, for very small detunings, CORPSE 
has more favorable scaling behavior. Figures 12a and 12b show the behavior of the shaped 



5 COMPOSITE PULSE SEQUENCES ON OTHER GROUPS 



40 



"w^--'' \ / ^ — / ''•■'.' 




-2-10 1 2 

Detuning error, eo 



Figure 13: a) Control function Ux(t) for various pulse sequences designed to compensate 
systematic detuning errors, b) Performance of several shaped pulses over a wide detuning 
range. 

pulse Si(7r/2) as trajectories on the Lie algebra and on the Bloch sphere respectively. 

5 Composite pulse sequences on other groups 

We have studied sequences which compensate imperfect single-qubit rotations, i.e., opera- 
tions which form a representation of the group SU(2). Another class of problems of practical 
and fundamental interest is the design of sequences for other Lie groups, such as the group of 
n-qubit operations SU(2") [50]. Several compensating sequences exist for multi-qubit gates 
[57, 83, 58]. In general, these sequences rely on operations which form an SU(2) subgroup 
to perform compensation in a way analogous to the one-qubit case and is the topic of this 
section. 

We note the design of compensating pulse sequences that do not rely on an SU(2) or 
SO (3) subgroup is a largely unexplored topic, and is an interesting subject for future study. 
For many cases, we can determine which errors can be compensated by examining the algebra 
[11], but it is unclear what the natural compensating sequences are for groups that are not 
equivalent to rotations in three- dimensions. 



5.1 Compensated two-qubit operations 

Thus far, we have shown how compensating sequences may be used to correct systematic 
errors in arbitrary one-qubit operations. Universal quantum computation also requires ac- 
curate two-qubit gates. The study of two-qubit compensating sequences is therefore of great 
fundamental and practical interest. In this section, we study two-qubit operations using Lie 
theoretic methods. The Cartan decomposition of the dynamical Lie algebra will be central 
to this approach. We begin by studying the properties and decompositions of the Lie group 
and its associated algebra. Then we consider systematic errors in two-qubit gates derived 
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from the Ising interaction, and show how the Cartan decomposition may be used to construct 
compensating sequences. 

Any two-qubit Hamiltonian (up to a global phase) may be written in the form, 

^W = EE^M^W^M^' (63) 

fj, V 

where the controls are represented in the product-operator basis [84], corresponding to the 
control Hamiltonians H^y = 2H^ ^ H^,, where Hi = ^1 and the indices /i, z/ run over 
Hfj^ e {Hi, ifx, Hy, Hz}. As a matter of convention, we exclude the term proportional to the 
identity Hu = ^1 (g) 1 since it generates an unimportant global phase and otherwise does 
not contribute to the dynamics; it is implied in the following equations that this term never 
appears. The product-operator representation is particularly convenient since each control 
Hamiltonian is orthogonal under the Hilbert- Schmidt inner product, {H^^, Hp^) = 5p_^p5y^„. 
The family of all possible solutions to a control equation, for example Eq. (1) with the 
Hamiltonian Eq. (63), forms a representation of the special unitary group SU(4). This 
group contains all possible single-qubit (local) operations that may be applied to each qubit 
among the pair as well as all two-qubit (nonlocal) operations. 

Cartan decomposition of two-qubit gates: Associated with the group SU(4) is the corre- 
sponding Lie algebra su(4) = span{— iif^^} excluding Hu (including Hu would make 
the group U(4)). Consequently su(4) is the algebra of four-dimensional traceless skew- 
Hermitian matrices. Observe that su(4) may be decomposed as su(4) = 1^ © m, where 

t = 5u(2) ® su(2) = span{-i/fxi, -i-^^yi) -i-f^zi, -i-f^ix, -i-^^iy, -i-f^iz} 
m = span{— iifxx, — i-f^xy, — i-f^xz, — i-f^yx, —^Hyy, — iiJyz, — iifzx, — i-f^zy, — i-f^zz}- 

It may be verified that [t, t] C t, [m, t] = m, and [m, m] C t. Therefore, the decomposition 
J^©m is a Cartan decomposition of the algebra su(4) (see section 3.3.3). Also, the subalgebra 
a = span{— iifxx, ~^Hyy, —iH^z} is a maximal abelian subalgebra and a C m. Therefore, a is a 
Cartan subalgebra. This admits the decomposition U = K2AK1 for any element U G SU(4). 
The operators 

ir,=exp|-i Y. «^?^Mi + «S^iJ (64) 

\ MG{x,y,z} J 

are in the subgroup e^ = SU(2)(8)SU(2) comprising all single-qubit operations for both qubits, 
whereas the operator 

A = exp -i ^ ap,pHf,p (65) 

\ /ie{x,y,z} J 

is in the abelian group e° = T^, isomorphic to the 3-torus, generated by the two-qubit inter- 
action terms ifxx, -f^yy, and H^^. In analogy with the Euler decomposition, the parameters 
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a^i, may be regarded as rotation angles. This construction is a KAK Cartan decomposition 
for the propagator U, and may be used as a framework for a pulse sequence to generate any 
arbitrary two-qubit operation. Furthermore, if the operators K2, A, and Ki may be im- 
plemented using a compensating pulse sequence, robust for a given error model, then their 
product will also be a compensating sequence. 

5.1.1 Operations based on the Ising interaction: 

We seek a method for implementing an arbitrary gate A G e" using a two-qubit coupling 
interaction. To be specific, we consider an NMR quantum computer of two heteronuclear 
spin qubits coupled by an Ising interaction. It is well known that propagators generated by 
the Ising interaction and single-qubit rotations are universal for SU(4) [85]; for completeness 
we explicitly show how arbitrary gates may be synthesized using the KAK form. Under the 
assumption of the qubits can be spectrally distinguished due to large differences in Larmor 
frequencies, the system Hamiltonian takes the form 

H{t)=l Yl u^i{t)H^,i + ui^{t)HiA +u,.{t)H,., (66) 

\M6{x,y,z} / 

where the controls u^iit) and ui^(t) are applied by the appropriate rf fields (see section 2.2) 
and Uzz{t) = 2nJ{t) is the strength of the spin-spin coupling interaction. In some cases, it is 
useful to manipulate the strength of the scalar coupling J{t), e.g., by using spin decoupling 
techniques [25] . We consider the simpler case of constant Ising couplings, and also assume 
that it is possible to apply hard pulses, where rf coupling amplitude greatly exceeds J, and 
the spin- spin coupling may be considered negligible. In terms of resource requirements for 
the quantum computer, this implies that single qubit operations are fairly quick, whereas 
two-qubit operations driven by the Ising coupling are much slower. 

Consider the propagators Ufj,^{a) = exp(— iaif^^) G e°. Let us define the one-qubit 
rotation operators Ri{9,(j)) = R{9,(t>) 1 and R2{9,(f>) = 1 (g) R{9,(t>). In the absence of 
applied rf fields, the system evolves according to t/zz(^zz), where 9zz = ^zzAt and At is the 
duration of the free-precession interval. Observe that for 

K^ = i?i(7r/2,0)i?2(7r/2,0)=exp(-i7r/2(if,i + /JiJ) 

Ky = i?i(7r/2,7r/2)i?2(vr/2,7r/2)=exp(-i7r/2(/7yi + i7iy)), (67) 

f4x(tt) = KyUzz{o:)K^ and Uyy{a) = KlXJzz{(y)Ky,. Since the group e° is abelian, then any 
arbitrary group element A specified by the decomposition angles {axx, ctyyjCtzz}, may be 
produced by the product A = U^{axx)Uyy{ayy)Uzzi(yzz) ■ Then any U = K2AK1 G SU(4) 
may be produced using single qubit rotations and the Ising interaction. 

Systematic errors in Ising control: In practice, systematic errors introduced by exper- 
imental imperfections prohibit the application of perfect two-qubit gates. In the case of 
gates produced by the Ising interaction, the errors may arise from several sources, such as 
experimental uncertainty in the strength of the coupling J . It is therefore desirable to design 
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sequences that implement accurate Ising gates over a range of coupling strengths. The error 
model in this case is similar to the case of amplitude errors; the control u^z is replaced by the 
imperfect analogue v^z = u^zi^ + ^j) where the parameter ej is proportional to the difference 
between the nominal (measured) and actual coupling strengths. For now we assume the 
remaining controls are error free. In correspondence, the perfect propagators Uzz{Ozz) are 
replaced by Kz(6'zz) = Uzz{Ozz{'^ + ej))- 

Jones was the first to study compensating pulse sequences for the Ising interaction 
[57, 55], and proposed sequences closely related to the Wimperis sequences studied in sec- 
tion 4.2. Here we demonstrate that these sequences are easily derived using Lie algebraic 
techniques. Observe that there are several subalgebras contained in su(4), for instance 
j = span{— iifxi, —iHyz, —iHzz} that are representations of 5u(2). Furthermore, if it is fea- 
sible to produce any imperfect propagator in the group J = e', then the compensating 
sequences discussed in section 4 may be reused for this system. 

This strategy can be implemented using accurate one-qubit rotations. For example, let 

n{e,(f)) = R\{(j), 0)Uzz{e)Ri{(j),0) = exp{-ie{cos(j)Hzz + sin (f)Hyz)). (68) 

Also, let us define the imperfect rotation A^( 61, 0) = i?|(0, O)Kz(6')i?i(0, 0) = 7^(6'(l + eJ), 0). 
The two qubit unitaries TZ{6, 4>) are isomorphic to single qubit rotations TZ{6, 4>) and Jones 
used this similarity to construct an alternative sequence we refer to as B2-J [57, 58], 

-Mb2-j(^,0) = Mi7r,(l))Mi27r,3^)Min,(P)Mie,0) 

= R\{(P,0)Uzz{7r{l + ej))R,{(P,0)R{{34>,0)Uzz{27i{l + ej))R,{34>,0) x 

i?I(0, 0)f/,,(7r(l + ej))i?i(0, 0) f/,,(0(l + ej)) 
= Uzzie) + Oie'j) (69) 

where again = arccos(— ^/47r). If the sequence B2-J is used in place of the simple Ising 
gate Uzz{6) then the first and second-order effects of the systematic error are eliminated. In 
this manner, any number of sequences designed for operations in SU(2) may be mapped into 
sequences that compensate Ising gates, including the higher-order Trotter-Suzuki sequences, 
which produce gates at an arbitrary level of accuracy. However, we note that since two-qubit 
gates occur so slowly, the practical utility of very long sequences is not clear, especially 
in systems where the two-qubit gate time is comparable to the qubit coherence lifetime. 
Substantial improvements in the minimum time requirements may be possible with shaped 
pulse sequences or using time-optimal control methods [86]. 

In this method, accurate one-qubit gates are used to transform inaccurate Ising gates 
into a representation of SU(2). Naturally, it is unimportant which qubit among the pair is 
rotated to perform this transformation, i.e., the subalgebra j' = span{— iiJix, —iHzj, —iHzz} 
would serve just as well. Given a control with a systematic error and a perfect rotation 
operator that transforms the control Hamiltonian H^ to an independent Hamiltonian H^,, it 
is possible to perform compensation if iiT^, and H^, generate a representation of su(2) [58]. 

Correcting simultaneous errors: Accurate single-qubit gates are required to compensate 
errors in the Ising coupling by transforming Ising gates into the larger dynamical Lie group. 



5 COMPOSITE PULSE SEQUENCES ON OTHER GROUPS 



44 



Amplitude error, €a 
-1.0 -0.5 0.0 0.5 1.0 




Figure 14: Infidelity of a) B2-J b) B2-W and c) B2-WJ. The target rotation is Ut = 
Uzz{t^,0)- The dashed contour corresponds to an infidehty of 0.01, while the remaining 
contours are plotted at 10% intervals. 

If accurate single-qubit operations are not available, then propagators of the form Eq. (68) 
can no longer be reliably prepared. However, if a compensating pulse sequence may be 
implemented in place of each imperfect single-qubit rotation, then the effect of this error may 
be reduced. This procedure was used in section 4.3.2 to produce a concatenated CORPSE 
and B2 sequence robust to simultaneous amplitude and detuning errors. A similar strategy 
can be employed to produce accurate Ising gates in the presence of simultaneous spin-coupling 
(two-qubit) and amplitude (one-qubit) systematic errors [58]. Consider the control system 
Eq. (66) under the influence of these two independent errors. The imperfect propagator 
takes the form 



V{u{t)) = U{u{t) + eA5u,{t) + ej5u2{t)), 



(70) 



where eA^Uiit) = e^Ux(i) + ^AUyit) is the amplitude error of the single qubit controls, 
and ejSu2(t) = eju^^it) is the error in the Ising coupling. The systematic error on the 
one-qubit controls e^<5iti(t) complicates the synthesis of compensated Ising gates, as unitary 
propagators of the form Eq. (68) can no longer be reliably prepared, i.e., Ai{6^ 0) = i?|(0(l-|- 
ca), O)[/,,(^(l + ej))i?i(0(l + e^), 0) = n{e{l + ej), 0) + C(e^). This difficulty may be avoided 
if we use a B2 sequence to correct e^i, before correcting ej using B2-J. Let t/B2(6',0) = 
f/r ® 1 + C'(e^) represent the propagator produced by a B2 sequence (see section 4.3.2) for 
the target rotation Ut = R{0, (p) on the first qubit. The sequence B2-WJ is 



M 



B2-WJ 



(^,0) 



M^2(0,O)f/. 



7r(l 
Mjf2(0,O)f/,,(7r(l 



ej))MB2(0, 0) iV42(30, 0)f/,,(27r(l 
ej))MB2(0,O)f/,,(^(l + ej)) 



ej))MB2(30,O) X 
(71) 



This sequence replaces imperfect i?i(0(l + eA),0) pulses with the compensated rotation 
produced by the B2 sequence. When e^ = 0, B2-WJ scales as C(ej), and when ej = the 
sequence scales as 0{e\). 
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5.1.2 Extension to SU(2'^): 

This demonstrated for two-qubits can be naturally extended to compensate operations on 
a network of n qubits with single qubit operations and Ising couplings. In this case the 
dynamical Lie algebra for this system is SU(2"). Khaneja has identified a recursive Cartan 
decomposition for this group that allows any U G SU(2") to be written as a product of 
single-qubit (local) and two-qubit (nonlocal) operations [50] . Given a pair of sequences for 
both single-qubit and two-qubit operations which compensate a particular error, it is always 
possible to produce a sequence for any U G SU(2") using this decomposition. In particular, 
imperfect Ising gates may be replaced with a B2-J sequence, while imperfect single-qubit 
gates may be corrected using the methods described in section 4. The surprising result is that 
one needs only a single accurate control to compensate an unlimited number of uncorrelated 
but systematic errors [58]. In practice this concatenating scheme is expensive for the whole 
system, but it points to a method for minimizing the amount of calibration required since a 
few good controls can compensate nearby errors. 

Computation on subspaces: Several interesting proposals involve the encoding of logical 
qubits on subspaces of a larger Hilbert space, which may offer certain advantages over other 
encoding schemes. In decoherence free subspace schemes, qubits are encoded on a subspace 
which is decoupled from environmental noise sources [87, 88, 89]. Also, several theoretical 
proposals involve the use of only two-qubit interactions to perform quantum computation 
[90]. In these schemes, the control algebra is chosen to be sufficiently large to allow universal 
computation on a subspace. 

An interesting question is whether gates applied to an encoded qubit may be corrected 
by using a compensating sequence on the codespace [58]. Recall that any gate in SU(2"') may 
be decomposed as a product of single-qubit (encoded) and two-qubit gates; it is sufficient to 
consider these cases individually. We may reuse the sequences described in section 4 if the 
controls for the encoded gates are distorted by a similar error model. For example, given a set 
of two controls with correlated systematic errors Vfj,(t) = {l + e)u^{t) and fjy(t) = {l + e)u^(t), 
it is possible to perform compensation if the control Hamiltonians H^ and H^, generate a 
representation of su(2). In Ref. [58], it is shown that universal subspace computation can be 
compensated if the two qubit Hamiltonians are of the XY model, H = H^^ + Hyy, but only 
single qubit operations can be compensated if the two qubit couplings are of the exchange 
type, H = H^ + Hyy + H^^. 

6 Conclusion and Perspectives 

Recent advances in quantum information and quantum control have revitalized interest in 
compensating composite pulse sequences. Specifically for the case of systematic control 
errors these techniques offer a particularly resource efficient method for quantum error re- 
duction. As quantum information processing experiments continue to grow in both size and 
complexity, these methods are expected to play an increasingly important role. 

In this review, we have presented a unified picture of compensating sequences based on 
control theoretic methods and a dynamic interaction picture. Our framework allows us to 
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view each order of the error as a path in the dynamical Lie algebra, highlighting the geometric 
features. Correction of the first two orders has a natural geometric interpretation: the path 
of errors must be closed and the signed area enclosed by the path must be zero. 

The geometric method helps illuminate the construction of arbitrarily accurate composite 
pulses. Currently arbitrarily accurate pulse sequences are of limited use because as the length 
of the sequence increases other noise sources become important. For most experiments, 
decoherence and random errors limit the fidelity of second-order compensation sequences. 
Our review of the Solovay-Kitaev pulses resulted in a modest improvement of the sequence 
time by implementing a new geometric construction. Shorter sequences may be achieved 
through numeric operation and continuous controls. 

Finally, we note that CORPSE and related pulse sequences are similar to dynamically 
corrected gates [4], in that both remove coupling to an external field while performing an 
operation. Combining the methods here with the developments in dynamically corrected 
gates and dynamic decoupling could lead to pulse sequences robust again environmental and 
control errors. These operations will lower the initial error and in the end limit the resources 
required for quantum error correction. This will ultimately determine the feasibility of 
performing large quantum chemistry calculations on a quantum computer [91, 92]. 
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